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1.0 SUMMARY 


This document describes the analysis and use of L219 (EQMOD), a digital computer 
program to modify matrices according to specific instructions. The particular field of 
application of the program is to modify the theoretical equations of motion and load 
equations generated by the DYLOFLEX programs Equations of Motion, L217 (EOM), 
and Load Equations, L218 (LOADS), respectively. 

The equations of motion and load equation coefficient matrices must be formulated 
outside of L219 (EQMOD) and read as program input data from magnetic files. These 
matrices can then be modified according to specific card input instructions, which 
allow the user to: 

• Use scalar multipliers. 

• Replace or increment individual matrix elements. 

• Add sensor equations to the equations of motion. 

• Add the definition of the active control system to the equations of motion. 

• Replace the theoretical rigid body and control stability derivatives in the equations 
of motion with those calculated by FLEXSTAB or other external means. 

• Transform the equations of motion and load equations from the inertia axis system 
to the body-fixed axis system. 

• Prepare matrix coefficients in a form useable in the Random Harmonic Analysis 
Program, L221. 

• Prepare matrix coefficients in a form useable in the Linear System Analysis 
Program, QR, to include: 

• Equations of motion with and without Wagner lift growth functions 

• Equations of motion and load equations combined for a time history solution 

The modified equations of motion and load equation matrix coefficients are saved on 
magnetic files for interfacing with solution programs L221 (TEV156) and QR. 

Limitations imposed on the matrices are: 

• The equations must be generated using the inertia axis system and a straight and 
level reference flight condition. 

• The vertical and lateral gusts are uncoupled (i.e., there are at most three rigid 
body degrees of freedom represented in the matrix coefficients for either the 
vertical or lateral gust analysis). 


2.0 INTRODUCTION 


The computer program L219 (EQMOD) may be used as either a standalone program or 
as a module of a program system called DYLOFLEX, which was developed for NASA 
under contract NAS1-13918 (ref. 1). Because of the DYLOFLEX contract requirements 
developed in reference 2, a program was needed to modify the equations of motion 
matrices generated by L217 (EOM) (ref. 3) and the load equation matrices generated by 
L218 (LOADS) (ref. 4). 

Modifications to the equations are needed when externally developed stability 
derivatives, such as experimental, or those developed in FLEXSTAB (ref. 5), are to be 
incorporated in the theoretical equations of motion. Other modifications such as varying 
the dynamic pressure and velocity, may be desired for parameter studies. These can be 
performed by matrix scalar modifications rather than regenerating the equations of 
motion and load equations with new aerodynamics. 
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3.0 SYMBOLS AND ABBREVIATIONS 


With the exception of section 6.3 (card input), all the items that appear in this 
document are listed below. 


ai.bj 

Ail 

b 

{C 3 } 

{ C 3? 

C D 

C L 

<7 

Pm 

Pn 

C y 

c* 

D 

D 


Dummy coefficients used for the mathematical description of a 
typical control system 

Wagner indicial lift growth function coefficients 

Time variant coefficients (functions of <j>i, 6\, cj) used in the 
development of the equations of motion. 

Reference span, length 

Matrix of generalized forcing function coefficients 
Matrix of load coefficients due to the excitation forces 
Aerodynamic drag coefficient 
Aerodynamic lift coefficient 
Aerodynamic rolling moment coefficient 
Aerodynamic pitching moment coefficient 
Aerodynamic yawing moment coefficient 
Aerodynamic side force coefficient 
Reference chord, length 
Aerodynamic drag force 

Generalized structural damping associated with elastic 
coordinates 


F X , Fy, F z 
F X'» FyS F z / 

W 

g 

h 


External forces defined in the body-fixed axis system. 

External forces defined in the interia axis system 

Matrix of streamwise distances from points first encountering 
the gust to the points encountering the gust later 

Gravity constant 

Acceleration at a particular sensor location 
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!xx> Iyy> ^zz> 
^xy> Ixz» lyz 

K 


L 

fi 

[Mi], [M 2 ], [M3] 


W], [m 2 ]. fM 3 ] 
[M4], [M 5 ] 

[M4], [m 5 ] 

M 

M 


M x , My, M z 

f , My I , f 

m 

n 

P. q. r 

q 

q 

s \\ r 

s 


Ut = V T 


Airplane inertias 

Generalized structural stiffness associated with the elastic 
coordinates 

Aerodynamic lift force 

Aerodynamic rolling moment 

Generalized structural stiffness, damping, and inertia matrix 
coefficients 

Load coefficient matrices (nonaerodynamic) of the generalized 
coordinate displacement rate, and acceleration, respectively 

Generalized aerodynamic stiffness and damping matrix 
coefficients 

Load coefficient matrices (aerodynamic) of the generalized 
coordinate displacements and rates 

Airplane mass 

Generalized mass associated with the elastic coordinates 
External moments defined in the body-fixed axis system 
External moments defined in the inertia axis system 
Aerodynamic pitching moment 
Aerodynamic yawing moment 

Rigid body roll, pitch, and yaw rates defined in the body-fixed 
axis system 

Generalized coordinate 

Dynamic pressure 

Wing reference area 

Laplace variable 

True velocity 
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u, V, w 

u e a > u e s 
W, 

x, y, z 
x',y',z' 

x cg» z cg 


Y 


Airplane velocities defined in the body-fixed axis system 

Perturbation values of antisymmetric and symmetric elastic 
coordinates 

Airplane vertical velocity in 1-g flight 

Linear displacement in the body-fixed axis system coordinates 

Linear displacements in the inertia axis system 

Coordinates of the airplane’s center of gravity in the body-fixed 
axis system 

Aerodynamic side force 


«g 


Gust angle of attack 


«1 


Trim angle of attack 


Exponential coefficients for the Wagner indicial lift growth 
function 


4 > 


(t) 


Wagner indicial lift growth function 


m 


Forcing function matrix when accounting for the gradual 
penetration of the gust 


[£] 


Matrix of load coefficients due to the excitation forces when 
accounting for the gradual penetration of the gust 


lJ>hJ 


<y,o',v 


8 c 

«Mt) 


Modal displacements at the sensor location 

Airplane angular displacements defined in the body-fixed axis 
system 

Airplane angular displacements defined in the inertia axis 
system 

The amount of rotation experienced by the control surface 

The amount of control surface rotating commanded by the 
stability augmentation system 

Kiissner indicial lift growth function 
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4.0 ENGINEERING AND MATHEMATICAL DESCRIPTION 

The equations of motion developed in L217 (ref. 3) are represented in the form: 

[Mi ] {q} + [M 2 ] {q} + [M 3 ] {q} + [M 4 ] {q} * + [M5 ] {q} * $ = {C3} a g * 0 d) 

where: 

M, C - Appropriate matrix coefficients 1 

q, q, q = Generalized coordinates and their time deratives, including SAS degrees of 
freedom 

«g = Gust angle 

3> = Wagner function (equal one for no lift growth) 

0 = Kussner function (equal one for no lift growth) 

* Indicates convolution 

With gust penetration, the excitation function, { C3 } dtg, of equation (1) is frequency 
dependent and is defined as: 


{C 3 } a g = [0] cos (Cl {%}) - i [01 sin (U {fjJ) 


( 2 ) 


where: 

Cl = o>/V r p, spatial frequency (rad/unit length) 

{fj^j = Matrix of streamwise distances from points first encountering gust to the 
points encountering the gust later 

[0] = Contribution of the lifting panels due to the gust force at designated 

gradual penetration load stations 

Relating these matrices in a physical sense, [Mj], [M2 ], and [M3] are usually 
associated with the generalized structural forces and the active control system 
definition. [M4] and [M 5] are usually associated with the generalized aerodynamic 
forces and {C 3} with the generalized excitation (gust) force. 


^ These matrix coefficients can be either constant or nonconstant (frequency dependent) 
coefficients. 
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The load equations developed in L218 (ref. 4) follow the same format as the equations of 
motion. 

{LOAD}= [M j ] {q}+ [M 2 ] {q}+ [M3] {q} + [M4] {q} * 

+ [M 5 ] {q} * $ + ( {C 3 } a g * 0) (3) 

where: 

[Mf], [M2], [M3] = Load matrix coefficients of the generalized coordinate 

displacement, rate, and acceleration, respectively 

[M4], [M5] = Load matrix coefficients of the generalized coordinate rate 

and acceleration convoluted with the Wagner function 

{C3} = Load matrix coefficient of the excitation function convoluted 

with the Kiissner function 

With gust penetration, the excitation function, {C^jag, of equation (3) is frequency 
dependent and is defined in a manner similar to equation (2): 

{C3 } = [ 0 ] cos (£2 {f^ }) - i [0] sin (£2 (f^ }) ^ 

where: 

[0] = contribution of the lifting panels due to gust forces at designated gradual 
penetration load stations to aircraft loads. 

Relating these matrices in a physical sense, [Mi], [M2], and [M3] are usually 
associated with the load resulting from structural response; [M4] and [M5] are usually 
associated with the load resulting from aerodynamic response forces; and { C3 } with load 
resulting from the gust excitation force. 

EQMOD offers an analyst the capability to alter the matrix coefficients of equations (1) 
to (4). Section 4.1 discusses the option to replace the appropriate theoretical 
aerodynamic terms in the M4, M5 and C3 matrices with airplane stability derivatives 
obtained from external sources such as flight or wind tunnel data. Section 4.2 details 
the option of incorporating the equations describing an airplane’s active control system 
into the basic equations of motion and load equations (all M, M, C and C matrices). The 
option to scalar multiply matrices or replace or increment individual matrix elements is 
presented in section 4.3. The option to transform the equations of motion of equation (1) 
into a form from which a stability analysis can be performed by solving for the 
eigenvalues of the equation is explained in section 4.4. Finally, section 4.5 discusses the 
changes of the matrix coefficients of equations (1) and (3) to transform from the inertial 
axes to the body-fixed axes. 
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4.1 STABILITY DERIVATIVES 


Developing the equations of motion for straight and level flight in the inertia axis 
system will result in some rigid body generalized coordinates acting at a vehicle 
reference point, representing rigid-body forward displacement (x), vertical displacement 
(z), and pitch (6) for symmetric flight conditions and lateral displacement (y), roll (</>), 
and yaw (tfj) for antisymmetric flight conditions. Some of the generalized coordinates 
may also represent rigid rotations about various control surface hinge lines. Embedded 
in the generalized aerodynamic and gust matrix coefficients of [M4], [M5], and {C3} 
are the total airplane theoretical aerodynamic forces and moments due to the airplane’s 
rigid body motions, control surface deflections, and gust angle of attack. For a 
symmetric analysis, these total airplane forces and moments correspond to the 
airplane’s lift, drag, and pitching moment defined in the inertial axis system. Similarly, 
for an antisymmetric analysis, these forces and moments correspond to the airplane’s 
side force, rolling, and yawing moments. 

These total airplane forces and moments in the inertial axis system can be related to 
appropriate airplane stability derivatives, usually defined in a particular stability axis 
system. This relationship is defined in detail in appendix A. A summary of the 
relationship between the stability derivatives and the appropriate total airplane force 
and moment elements of [M4], [M5], and [C3] are shown in figures 1, 2 and 3. If these 
derivatives are available from wind tunnel results, flight test results, or any other 
source, of if they are calculated in FLEXSTAB (ref. 5), they can be used to calculate the 
appropriate matrix coefficients and used in preference to the theoretical coefficients 
calculated in L217 (EOM). 

Since elastic modal degrees of freedom (elastic generalized coordinates) are included in 
the dynamic analysis, the stability derivatives used in the equations of motion should 
be only rigid body stability derivatives. The aeroelastic effects that are represented by 
elastic increments to the rigid body stability derivatives are reflected in the equations 
of motion through the elastic modal representation. However, since DYLOFLEX in 
general does not consider panel aerodynamics that are not perpendicular to lifting 
surfaces, no aeroelastic effects are represented in the forward generalized coordinate 
displacement. Consequently, the symmetric representations in figures 1 and 3 use both 
the rigid stability derivatives and the elastic increment to the rigid stability derivative 
in calculating the generalized matrix coefficients for the forward (x) generalized 
coordinate to obtain aeroelastic effects for that degree of freedom in the dynamic 
analysis. 

In addition, the development of the expressions in figures 1 through 3 assumes that the 
stability derivatives and the rigid body motions are defined about the same reference 
point. If the reference points are at different locations, the corrections that must be 
made to the stability derivatives in figures 1 to 3 are shown in figure 4. 



Aerodynamic stiffness matrix [MJ 




-C O 
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Figure 1. — Formulation of The Rigid-Body Symmetric Generalized Aerodynamic Stiffness and 
Damping Matrix Elements Using Stability Derivatives 




Aerodynamic stiffness matrix [M a ] 



Figure 2. — Formulation of the Rigid-Body Antisymmetric Generalized Aerodynamic Stiffness and 
Damping Matrix Elements Using Stability Derivatives 





















(a) Symmetric 


Gust forcing column {Cgl 



This is the z^ql of the Mg matrix with the sign changed. (Note: Due to the sign convention adopted in EOM, 
only the symmetric case requires a sign change:) 


(b) Antisymmetric 

Gust forcing column IC-jl 



Figure 3. - Formulation of the Rigid-Body Gust Excitation Matrix Elements Using Stability Derivatives 
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(a) Symmetric 


Aero ref = 
Ref 


Cl V^ Cd 0r 

REF 


Ax A z 

C m = C m +— C, C n 

m or R m a R c L a R c D a 

REF 


Ax 


Az 


C mA C mA + C LA ‘ - C DA 
q R q R c q R c q R 

REF 


Ax Az 

C m “ G__A + — Q | ^ * - Cr> A 

m O' R m a R c L a R c D <* R 

REF 


Ax 


Az 


C_» = ^ Ci ~ Cr 


m r m 


c L « 5 0 c u 5 , 


REF 


Where Ax “ < X REF " x AERO REF* 
Az = (z REF ~ z AERO REF) 


+z, up 

^ — ► +x, aft 

Reference point about which the aerodynamic derivatives are calculated. 
Reference point about which the rigid-body motion is defined in the analysis. 


Figure 4. — Stability Derivative Corrections for a Different Aero Reference Point 



(b) Antisymmetric 


X 

- Co - 

Az 

b 

Cy 0 R 



REF 

Az 

TT 



X 

" CnA - 

X 



REF 

Az 



X 

= " 

P R 

X 

and 

r R 

REF 

= C 9a ~ 
^R 

Az 

b 

X 


% 

6 r 

REF 

= Co - 
5 R 

Az 

b 

X' 



REF 


s 

P R 

= c n + 
% 

Ax 

b 

X 


REF 

Ax 

b 


C,,A 

% 

= C n A + 
% 

X 

% 

REF 

= C "£ + 
P R 

Ax 

b 

X 

X 

REF 

= * 

r R 

Ax 

b 



REF 

Ax 

b 


X 


Cy *R 


REF 




Figure 4 , — (Concluded) 
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4.2 ACTIVE CONTROL SYSTEM DEFINITION AND SENSOR EQUATIONS 

If the effects of an active control system are to be represented in the equations of motion 
of an airplane (eq. 1), it is necessary to define the active control system as a number of 
linear second order or less differential equations. As an example of this procedure, a set 
of linear differential equations are developed using the sample control system in 
figure 5. This control system is only an example used for the purpose of illustration. The 
user is free to use any type of control system as long as its mathematical description 
can be reduced to a set of linear second order or less differential equations. 



Figure 5. — Sample Control System 


where: 


qj to q n 
Qn 
Sc' 
h 

Qn+l to 
Qn+3 


= Vehicle degrees of freedom (generalized coordinates) 

= Control surface rotation (S c ) 

= Amount of control surface rotation commanded by the control system 
= Acceleration at a particular sensor location 

= Dummy degrees of freedom 


Throughout this section, the control system equations will be developed in the Laplace 
domain (s-plane); therefore, all generalized coordinates and control surface rotations 
will be functions of the Laplace variable, s. 


In the example of figure 5, the equations of motion were derived using n degrees of 
freedom with the n^ generalized coordinate being the actual control surface rotation, 
8 C . However, this is not a requirement of the program. The control surface rotation may 
occupy any position in the generalized coordinate array. It is also assumed, for the sake 
of simplicity, that the amount of rotation experienced by the control surface, 8 C , will be 
equal to the amount of rotation commanded by the control system, 5 C '; that is, the 
control system is perfect. This is usually not the case. In reality the physical properties 
of a control system (e.g., the deflection of backup structure, maximum actuator force 
available, etc.) coupled with the aerodynamic forces on the control surface will result in 
a difference between the amount of control rotation commanded and the actual amount 
experienced. In such instances, the user must supply the appropriate equations that will 
describe this difference. 
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In the following equations, qi to q n represent the airplane’s elastic and rigid body 
degrees of freedom. The quantity h will be an acceleration sensed on the vehicle and 
used as feedback input by the control system. In general, the quantity h may be 
displacement, velocity, or acceleration. q n +i to qn+3 are dummy degrees of freedom 
used to keep the control system equations to a set of second order or less differential 
equations. 

Working through the block diagram shown in figure 5, the set of differential equations 
describing the sample control system is derived as follows: 

a) Based on the assumption of a perfect control system 

5c = 5c' 

or 

5c - 5c' = 0 (5) 


b) From the outer loop shown in figure 5 

_ a 

Qn+1 ~ s + b 


or 


(s + b) q n+1 - ah = 0 


( 6 ) 


c) The relationship between the acceleration, h, and the generalized coordinates 
excluding the dummy coordinates is 


= h 


U>h J - 

% . ' 

Ql 

• 

• 

• 

’ = S 2 L0 h J " 

V 

• 

• 


Qn 



or 

r 




s 2 L0 h J 


Ql 

. l-h = 0 

Qn 


(7) 


The mode shapes L^J are the modal deflections at the particular sensor location 
and they can be obtained from the program LOADS(L218). Equation 7 may be 
considered as the sensor equation. 


d) From the summation point 


Qn+2 Qn+1 Qn+3 
or 

Qn+2 “ Qn+1 + Qn+3 = 0 


( 8 ) 
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e) The forward path gives 


6c' 


+ ds + e 


9n+2 


or 


(s^ + ds + e) 5c' - cq n+ 2 = 0 
f) Finally, from the inner feedback loop 


f 

9n+3 s + g 


6c' 


( 9 ) 


or 

(s + g) q n +3 = 0 (10) 

The coefficients for equations (5) through (10) relating the addition degrees of freedom 
(q n+ l, 9n+2> 9n-l-3.Sc' and h) to the unaugmented vehicle degrees of freedom (qj to q n ) 
must be added to the unaugmented equations of motion matrix coefficients. 


Taking’ the Laplace transform of the equations of motion (eq. 1), the unaugmented 
equations of motion matrix coefficients of [M^], [M2], and [M3] represented by square 
matrices of size (n X n), can be written as coefficients of s. 



(ID 


The equations of motion are now expanded by five degrees of freedom with the 
coefficients of equations (5) through (10) placed in the appropriate matrix locations as 
shown in the following: 


Mi 

J UA 

1 

0 

-1 


f ^ 

<11 

• 

• 

9 n -l 

5c 


M 2UA 

0 


91 

• 

• 

9 n -l 

5c 


b 0 0 -a 0 

< 

9n+I 

• + s 


1 0 0 0 0 


9 n +l 


-11 100 


9 n +2 



00000 


9 n +2 

0 

0 -c 0 0 e 


9 n +3 


0 

OOOOd 


9 n +3 


0 0 0 -1 0 


h 



0 0 0 0 0 


h 


0 0 g 0 -f 


Sc 

L j 


_ 

0 0 10 0 


5 c' 

V. y 
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Mo 

3 UA 

0 


^1 

• 

• 

• 

%-l 

5 c 


0 0 0 0 0 

< 

%+l 

0 

0 0 0 0 0 


%+2 


0 0 0 0 1 


%+3 

L <£ h J 

0 0 0 0 0 


h 

0 

0 0 0 0 0 

J 


5 c' 

l 


( 12 ) 


Each added row in equation 12 corresponds to one of the equations 5 through 10 . 


The inclusion of the control system equations, as illustrated in this example, is 
accomplished in EQMOD by the user specifying on cards the coefficients and their row 
and column locations in the augmented matrices. The aerodynamic related matrices, 
[M4], [M5] and either {C3} or [<£>], are automatically augmented by rows and columns 
of zeroes. No aerodynamic forces and moments result from adding the dummy degrees of 
freedom. 

It is important to note that the sensor equation for this sample control system used 
acceleration as a feedback quantity. Therefore, the mode shape matrix, </> n was placed 
in the augmented M3 matrix. In the case of displacement or velocity being the feedback 
variable, the modal data would have been placed in the augmented M] or M2, 
respectively. These sensor coefficients (mode shapes) may be either calculated in L 219 
(LOADS) (ref. 4 ) and read directly into L 219 (EQMOD), or they can be input manually 
into the program by cards. 

Also with regard to the sensor equation ( 7 ), the variable h is automatically included by 
EQMOD in the augmented generalized coordinate array. In defining the sensor data, 
the user need only specify the matrix (Mj, M2, M3) and the row at the matrix where the 
mode shapes are to be placed and EQMOD automatically places a -1 in the groper 
column of Mj. In the sample problem, the mode shapes were placed in the (n + 4 ) u row 
of thegjugmented M3 matrix (eq. 12), therefore the -1 was placed in the (n+ 4 ) ul row and 
(n+ 4 ) column of the augmented M^ matrix. 

This example used only one sensor h. EQMOD can accept more than one sensor and 
more than one type of sensor. However, the method used for more sensors is identical to 
that used in this example. 
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4.3 MATRIX MODIFICATION BY SCALAR MULTIPLICATION, 

REPLACEMENT OR INCREMENTATION OF MATRIX ELEMENTS. 

EQMOD offers the capability for modifying all the coefficient matrices in the equations 
of motion (1) and (2) load equations (3) and (4). The matrix coefficients may be modified 
by multiplying the entire matrix by scalar factors or by replacing or incrementing 
individual elements within the matrices. The following is an example of matrix 
modification by scalar multiplication. 

In the generation of the [M4] and [M4 matrices, the dynamic pressure is embedded in 
the coefficients and thus [M4] and [M4] are functions of the dynamic pressure, q. 

[M 4 ] , [M 4 3 =/(q) (13) 

Similarly, the [M 5 ], [M 5 ], {C 3 }, or [<?] and {C 3 } or [$] matrix coefficients are 
functions of q and the reciprocal of the freestream velocity, Vt. 

[M 5 ], [M 5 ], a {c 3 } , {c 3 } = /( q ‘ ) (14) 

If a parameter study is desired varying only the dynamic pressure and freestream 
velocity, but keeping everything else constant including the Mach number, the [M4] 
and [M4] matrices can be multiplied by a scalar which is the ratio of the new dynamic 
pressure to the original dynamic pressure. 

r w i _ q new rx . , 

[M^new q [ M4I 


[M4]new=^ IM 4 ] <1£ 

[M 5 ], [M 5 ], {c 3 } or [?], and {C 3 } or [<£] matrices can be multiplied by a scalar 
consisting of this ratio times the ratio of at the original velocity to the new velocity 


[^5] n ew 


Vj 

Vj 


new 


q new 

q 


[m 5 ] 


(16) 


and so on. 

In performing the matrix modifications in EQMOD, the values of the scalar multipliers, 
the replacement elements, and the element increments must be calculated by the user 
and input via cards. EQMOD dos not calculate those values internally. 
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4.4 FORMATION OF EQUATION OF MOTION CHARACTERISTIC EQUATION 
WITH WAGNER INDICIAL LIFT GROWTH FUNCTION 


A method of obtaining the stability of a system is to calculate the roots of the system’s 
characteristic equation. The Linear System Analysis Program (QR) (ref. 6) has the 
capability of calculating the characteristic equation of a system from the Laplace 
transform of the equations of motion (1) and determining the roots of that equation. It is 
the function of EQMOD to formulate a set of matrix coefficients which represent the 
Laplace transform of equation (1). The following is the theory to form this set of matrix 
coefficients that includes indicial lift growth functions applied to the equations of 
motion generalized aerodynamic coefficient matrices. 


In the time domain, the Wagner indicial lift growth function may be approximated as 

-Of 1 1 -/3 1 1 

$(t) = 1 - aje 1 - b je 1 (171 

The Laplace transform of equation (17) is 

1 a l b l 

^(s) = — t~~ - c , V (is: 

In equation (1), the Wagner indicial lift growth function is convoluted with q and q. 
Using the relationship for the Laplace transform of the convolution integral (Duhamel’s 
formula): 

y [ / l f(t - x> g (X) dx] = <e [f(t) * g(t)j =y[f(t)i y[g(t» <19 

o 

the left side of equation (1) in the Laplace domain becomes: 

[[Mj] + s[M 2 ] + s 2 [M 3 ] + s[M 4 ] -1(1 - .iiL.) 


+ * 2 i M si f (1 - rr^ t ■ 0 < 20) 

Expanding and collecting terms gives: 

[s 4 [M 3 ] + s 3 [[M 2 ] + [M 5 ] + (aj + j3 1 )[M 3 ] - (a, + bj)[M 5 l] 

+ s 2 [[Mj ] + [M 4 ] + (aj + j3 1 )[[M 2 l + [M 5 ]] + aj /3 j[M 3 ] 

- (a + b)[M 4 ] - (a j/3j + b 1 a 1 )[M 5 ]] 

+ s^(aj +0i)jjMj] + [M 4 lj + a l^l[fM 2 ] + [M 5 ]] 

-(a j/Jj +b 1 of 1 )[M 4 ]] +« 1 i 8 i[[M 1 l + [M 4 ]]j {^(q)} = 0 ( 21 ) 
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If indicial lift growth is not included, that is: 

* (t) = 1 
or 

$ (s) =y (22) 

and consequently aj and bj from equation (18) are zero, equation (21) will simplify into 
only a second power of s equation. 

The characteristic equation of the system is obtained by taking the determinant of 
equation (21) and setting it equal to zero. It should be noted that due to the inclusion of 
the Wagner function, the order of the characteristic equation is increased and additional 
roots will be calculated. These additional roots do not represent additional modes of the 
system. The number of additional roots is a function of the order of the set of 
polynomial equations shown in equation (21) and of the number of degrees of freedom in 
the system. The order of the set of polynomials is a function of the number of terms 
used in the approximation of the Wagner function. 

The function of EQMOD is to form the matrix coefficients of s shown in equation (21). 


4.5 TRANSFORMATION FROM INERTIAL AXES TO BODY-FIXED AXES 


The matrices for equations (1) and (3) must be generated in the inertial axis system 
for straight and level flight. If the user prefers to work in the body axis system, a 
transformation from the inertia to body-fixed axis system may be desirable and can be 
accomplished. 

Presented in this section is a summary of the matrix changes that are made by EQMOD 
to the coefficients of equations (1) and (3) in order to convert from inertia axes to 
body-fixed axes. A full theoretical development of the transformation is given in 
appendix A. 

Basically, the effect of the transformation is to redefine the generalized coorinates that 
describe the rigid body motions of the aircraft in equations (1) and (3). All other 
coordinates (elastic and control deflections) are not affected. In the inertia axes, which 
are fixed in space, the motion of the aircraft relative to these axes is described by the 
velocity components in the direction of the inertia axes. In the body-fixed axes, however, 
the motion is described by the velocity relative to the fixed inertia axes but in the 
direction of the moving axes. The relationships between the velocity components in the 
inertia axes and the body-fixed axes are shown in table 1. 


U] and Wj are the reference (in this case lg) values of velocity defined in the reference 
axis system which is fixed. Uj will be referred to as the airplane forward velocity, Wp. 
For small angles, the reference (lg) angle of attack can be defined as 


<*} ~ tan <*1 



(23) 
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Table 1. — Relationships Between Velocity Components in the Inertia Axes 
and the Body-Fixed Axes 


Analysis 

Inertia axes 

Body-fixed axes 


* / 
X 

u + W-,0' 

Symmetric 

z' 

w - U ^0 ' 


Q ' 

q 


y' 

v + U.,0'- W.,0' 

Antisymmetric 

0' 

p 


i' 

r 




The effects of the transformation on the meaning of the generalized coordinates and the 
changes in the various matrices will be examined for the symmetric analysis first and 
the antisymmetric analysis second. 


Symmetric Analysis 

In the symmetric analysis, the generalized coordinates can be interpreted as: 
In the inertia axis system 



r •> 

X' 


. X 

x' 


fx'l 

{q} = - 

z' 

6 ' 

• {<]} = • 

z' 

d’ 

' {q} -< 

z" 

6 ' 


% 


^e 


% 


L 5 cJ 


5 C 

L I* J 


5 c 

k. ^ 


and in the body-fixed axis system 


r ' 

X 


r 

u 



y 

> {q} - ’ 

w 


w 

e 

q 

% 

■ {q} = • 

q 

% 

5 c 

V ^ J 


5 c 

L ^ J 


K 

i c j 


(24) 


(25) 


Note the elastic coordinates, q e , and the control surface coordinates, 5 C , remain 
unchanged from one system to the next. 


To transform equations (1) and (3) to body-fixed axes, the rigid body generalized 
coordinate velocity and acceleration matrices in the inertia axes are replaced with the 
expressions given in table 1 and with derivatives of these expressions. The resulting 
terms are then regrouped into coefficients of generalized coordinate displacements, 
velocity, and acceleration. In doing so, the transformation from inertia to body-fixed 
axis system requires the following changes to the coefficient matrices [Mj], [M2], 
[Mi], [M|], [M2], and [Mj], 


The only column changed in these matrices is the 6: 



^col 

to 


M, 

'•col 

- v T ( 

m 2 

z col 

* a l M 2 ) 

x coi ; 

Mq 
2 « col 

- V t( 

M3 

z col 

* a l M 3 ) 

x col 1 

M 4 S 

^col 

- v T ( 

M 5 z 

-«1 m 5 ) 


z col 

x col ^ 


in[Mj] : 
in[M2] : 
in[M 4 ] : 

; w col ' ^col A col ' 

M n = the a column of the original (inertia axis) matrix of the n^ matrix 
a col 


(26) 


[M,], [M2], and [M4] are changed in the same manner. The transformation does not 
affect the [M3], [M3], [M5], or [M5] matrices. 
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Antisymmetric Analysis 

Similarly, in the anti-symmetric analysis, the generalized coordinates can be 
interpreted as: 

In the inertia axis system: 


{q} - ' 

(VI 

r 

{q} = - 

[y'l 

4>' 

f 

• ' {q} = < 

f y' 

r 




% 

L*cJ 


% 

5 c 


And in the body-fixed axis system: 



fy ] 
<t> 


r v ' 
p 


r • 

V 

P 

{q} - ■ 

* 

% 

5 c 

L ^ J 

{q} = ' 

r 

-{q} = • 

r 

% 

s c 


Performing a similar substitution as with the symmetric analysis — but with the 
antisymmetric variables in table 1 — and regrouping, the following changes to the 
coefficient matrices are required: 


*col to 'j'col to 


in [Mj] 

Mi - V T atMn 

^ col y col 

Mi , + V T M 9 

^coi y C oi 

in [M 2 ]: 

M 2 - V T a l M 3 v . 

coi y coi 

M 9 + VtMo 

^coi ycoi 

in [M 4 ] : 

- Vt QIiMc 

4* col V col 

%coi +VTM5 y C oi 


[Mi], [M2], and [M4] are changed in the same manner. Again, [M3], [M3], [M5], and 
[M 5 ] are not affected by the transformation. EQMOD performs all matrix 
manipulations needed to complete the transformation. The user needs only define the 
column locations of the rigid body motions, the Ig angle of attack and the airplane’s 
forward velocity, Vq\ 
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5.0 PROGRAM STRUCTURE AND DESCRIPTION 


L219 (EQMOD) has been constructed as an overlay system. Figure 6 shows the overlay 
structure and the data input to and output from each overlay. The overlays are: 


Main overlay (L219,0,0) L219vc 

Primary overlay (L219,l,0) RDCRDS 

Secondary overlay (L219,l,l) RDEOM 

Secondary overlay (L219,l,2) RDLOD 

Secondary overlay (L219,l,3) RDQR 

Primary overlay (L219,2,0) EOMMOD 

Primary overlay (L2 19,3,0) LODMOD 

Primary overlay (L219,4,0) QRMOD 


The main overlay, L219vc ("v” and "c” are version and correction identifiers), is a small ' 
program that calls into execution the primary overlays required to perform the^ 
operations requested by the user via card input data. 

The first overlay called into execution by L219 is always the 1,0 primary overlay 
RDCRDS. RDCRDS reads and interprets all card input data, determines the execution 
options, and writes the edited input data onto the random scratch file SCRAND. 
RDCRDS calls three secondary overlays to process special sections of the input data: 

CL21 9, 1,1) RDEOM reads instructions from cards directing the modification of 

equations of motion matrices. FLEXSTAB (with the DYLOFLEX 
modification, ref. 7) stability derivatives are read from the file SDSSTP, 
if required. Also, sensor equations will be read from the file LODTP2, if 
requested. 

(L219,l,2) RDLOD reads instructions from cards directing the modification of load 
equation matrices. 

(L219,l,3) RDQR reads instructions from cards directing the preparation of 
matrices for the program QR. 

The remaining primary overlays (2,0, 3,0, and 4,0) will be executed only if requested by 
the user. All primary overlays read input instructions from SCRAND. 
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Overlay 2,0 (EOMMOD) modifies the equations of motion matrices read from the file 
EOMTAP. The modifications are made in the following order (all are optional): 

• Calculate stability derivatives and store over the proper matrix elements. 

• Scalar-multiply matrices. 

• Add sensor equations. 

• Replace and increment elements. 

• Perform the inertia to body*axis transformation. 

The resulting matrices are written onto the file EQEOM specified by card input in 
Overlay 1,0 for use in the Random Harmonic Analysis program, L221 (ref. 8). 

Overlay 3,0 (LODMOD) modifies the load equation matrices read from the file 
LODTAP. The changes are made in the following order. 

• Scalar-multiply matrices. 

• Replace and increment elements. 

• Perform inertia to body-axis transformation. 

The resulting matrices are written onto the file EQLOD specified by card input in 
Overlay 1,0 for use in the Random Harmonic Analysis program, L221, (ref. 8). 

Overlay 4,0 (QRMOD) prepares matrices for the program QR. The necessary equations 
of motion and load equation matrices are used to form the matrices for rooting and for a 
time history solution. The matrices are written onto the file QRTAP as specified on card 
input in Overlay 1,0. 

For a more complete description of the L219 (EQMOD) program structure see volume II 
of this document. 
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6.0 COMPUTER PROGRAM USAGE 


The program was designed for use on the CDC 6600. The machine requirements to 
execute L219 (EQMOD) are: 

To read control cards and card input data. 

To print standard output information, optional intermediate results and 
diagnostic messages. 

All magnetic files not specifically defined as magnetic tapes are 
assumed to be disk files used for input, output, and temporary file 
storage. 

For "permanent” storage of data; magnetic files are copied to and from 
magnetic tapes with control cards before and after program execution. 

The program L219 (EQMOD) is written in FORTRAN and may be complied with either 
the RUN or FTN compiler. L219 may be executed on either the KRONOS 2.1 or NOS 
operating system. 


Card reader 
Printer 

Disk storage 
Tape drive 


6.1 CONTROL CARDS 

The following list is a typical set of control cards used to execute L219 (EQMOD) using 
the absolute binaries from the program’s master tape. 

Job card 

Account card 


REQUEST(MASTER,F=I,LB=KL,VSN=66XXXX) 
RE WIND( MASTER) 

SKIPF(MASTER) 

COPYBF(MASTER,L219) 

RETURN ( M ASTE R) 


Retrieve the 
< program from its 
Master tape 
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Prepare optional 
input data files 


L219, 


EXIT. 

DMP(0, field length) 
— End-of-record 



| Excecute L219 (EQMOD) 

( Save optional output 
data files 


6.2 RESOURCE ESTIMATES 

The computer resources used (core requirements, tapes, printed output, time, etc.) are a 
function of the problem size and program options used. 

FIELD LENGTH 


The field length required by L219 (EQMOD) is dependent upon the problem size and the 
program module(s) used. Core must be requested based on the largest amount of core 
needed for any one module to be run, that is: 


Field length 


Maximum required for 


f RDCRDS 
EOMMOD 
LODMOD 
^ QRMOD 


For each module, the core requirements are determined from the following formulas 


RDCRDS 110,000g + NDOF*NDOF 

EOMMOD 71,000(8 + 3*NDOF*NDOF 

LODMOD 66,000(8 + 2*NODF*NDOF + NDOF*NLD 

QRMOD 67,000(8 + 2*NDOF*NDOF 
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where: 


NDOF = number degrees of freedom 


NLD = number of loads 
Time Estimate 

The time estimate is dependent upon the problem size. However, the average time to 
run most average jobs should be less than 25 seconds. 

Printed Output 

The maximum number of lines of printed output has been limited to 40 000, which 
should be enough for any L219 (EQMOD) program execution. The average line count is 
about 1000 lines. If output line limit is exceeded, use the following control card to 
execute the program: 

L219(PL = limit) 

where "limit” is the approximate number of lines required to execute this program. 

6.3 CARD INPUT DATA 

A detailed description of the card input data needed to execute EQMOD is given in 
sections 6.3.1 through 6.3.4. A summary of the card input data is given in section 6.3.5. 
The summary is a quick reference for the necessary card input and is included for use 
only after familiarity with the program has been obtained. 

The task(s) performed by L219 (EQMOD) are broken into three subtasks, each with its 
own section of code known as a primary overlay. The entire set of primary overlays is 
driven by a small program (main overlay) named L219vc. 

L219vc reads program directive cards to: 

• Assure that the data being read is intended for L219 (EQMOD). 

• Determine which section of code (primary overlay) of L219 is to be executed next. 

• Determine what data and results are to be printed. 

• Determine input and output magnetic File names. 

• Determine total problem size. 

The order in which the input cards are read is shown in figure 7. 
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<>-* 


f 3 

\ PRINT INPUT 


I 

Q*+ 


*».3 


PRINT OUTPUT 


(^SYMMETRIC 

IANTI- 

(SYMMETRIC 



Figure 7. — Flow of Card Input Data 


Comment cards may appear 
anywhere in the input 
stream except where data 
follows a keyword card 
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LOADS EQUATIONS (Repeated for each file) 










Format of Card Input Data 

All card data is read in fixed fields, specific columns of the cards. On the pages that 
follow, the required card columns are defined next to each keyword or variable. The 
following conventions, which are used throughout the program, should be noted. 

• All floating point variables are read with format E10.0. 

• All integer variables are read with the format 15. 

• All hollerith variables (keywords, etc.) are read with the format A10; however, 
when the program is trying to recognize keywords, it checks only the first four 
characters. 

All data fields end on a card column that is a multiple of five. 

6.3.1 GENERAL OPTIONS 

The first card read by L219 (EQMOD) must be $EQMOD, card set 1.0. It indicates that 
card data for L219 (EQMOD) follows. 

After Card Set 1.0 the program continues to read data cards and checks the first four 
characters for keywords. The keywords introduce the remaining card sets. Card sets 2.0 
through 5.0 define the problem size, file names, and options to be used throughout the 
program execution. 


Card Set 1.0 - Equation Modifier (EQMOD, L219) for 

Equations of Motion and Load Equations 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

$EQM Od 

A4 , 6 X 

This card must be the first card read by the Equation 
Modifier (EQMOD) program. The $EQMOD card indicates 
that the data following is for the Equation Modifier 
Program. 


Card Set 2.0 - Case Labeling Information 


Card 2.1 - Title Card 

The title card after the $EQMOD card will be stored in core, up to four title cards, for 
page headings on printouts. 
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Card 2.2 - Comment Card (Optional) 


Comment cards may appear anywhere in the input data stream except where data 
follows a keyword card. 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-2 


A2 

Keyword for comment card. 

NOTE; A blank in column 2 must follow the C in 
column 1. 

3- BO 

COMMENT 

A8,7A10 

Comments will appear in the printed output as they 
are read. It is not treated as data. 


Card Set 3.0 - Problem Size 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

SIZE 

A4 , 6X 

Keyword indicating the problem size. 

11-15 

NDOF 

15 

Total number of output degrees of freedom (no default) 
(<_ 100) 

16-20 

NPAN 

15 

Total number of output panels {number of gust zones) 

NOTE: If the equations of motion are generated external 

to DYLOFLEX set 

NPAN = 0 if no gradual penetration is being used 

and is real. 

(Default: NPAN extracted from EOM tape) (<^ 50) 

21-25 

nfreqm 

15 

Number of frequencies at which unsteady aerodynamics are 
defined. (Default: NFREQM extracted from EOM tape) 

<i 20) 
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Card Set 4.0 - Output Options 
Card 4.1 - Output Tapes 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

OUTPUt 

A4,6X 

Keyword indicating output file names and file positions. 

11-20 

IUTEOM 

A7, 3X 

File name where equations of motion matrices are written. 
(Default: IUTEOM = EQEOM) 

21-25 

IPLEOM 

15 

File position number where equationsof motion matrices 
are written. 

(Default: IFLEOM = 1) 

26-30 

dummy 

5X 

Blanks . 

31-40 

IUTLOD 

A7,3X 

File name where load equations matrices are written. 
(Default: IUTLOD = EQLOD) 

41-45 

IFLLOD 

15 

File position number where load equationsmatrices are 
written. 

(Default: IFLLOD = 1) 


Card 4.2 - Print Input Matrices (Optional) 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

PRINt 

A4 , 6 X 

Keyword introducing the print option. 

11-20 

INPUt 

A4,6X 

Keyword indicating the input matrices to be printed. 

21-30 

MATRIX 

10X 

Descriptive word (optional) 

31-40 

OPTION 

A10 

Print options: 

Keyword Matrices Printed 

ALL All input matrices printed 

NONE No input matrix printed 

FREQUENCY Only input matrices of Ith frequency 

(ITHF) printed. 

Default: (OPTION = NONE) 

41-45 

ITHF 

15 

(Required only if OPTION ■ FREQUENCY) 

Ith frequency input matrices to be printed. 
(Default: ITHF = 1) 
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Card 4.3 - Print Output Matrices (Optional) 


COLS. 

1-10 

11-20 

21-30 

31-40 


41-45 


KEYWORD/ 

VARIABLE 


PRINt 

OUTPut 

MATRIX 

OPTION 


ITHF 


FORMAT 


A4 , 6X 
A4,6X 
10X 
A10 


15 


DESCRIPTION 


Keyword introducing the print option. 

Keyword indicating the output matrices to be printed. 
Descriptive word (optional) 

Print options: 


Keyword 

ALL 

NONE 

FREQUENCY 


Matrices Printed 


All output matrices printed 
No output matrices printed 
Only output matrices of Ith frequency 
(ITHF) printed. 

CHANGED Only those matrices that have been 

changed are printed. 

(Default: OPTION = CHANGED) 

(Required only if OPTION = FREQUENCY) 

Ith frequency output matrices to be printed. 
(Default: ITHF = 1) 


Card Set 5.0 - Symmetric or Antisymmetric Analysis (Optional) 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

fSYMMetric 1 

A4 , 6X 

Keyword SYMMETRIC indicates a symmetric analysis for 
body axis and derivatives. 

Keyword ANTISYMMETRIC indicates an anti-symmetric analysis 
for body axis and derivatives. 

(Default: SYMMETRIC) 

[ ANTI symmetric | 
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6.3.2 INSTRUCTIONS TO MODIFY EOM MATRICES 

Omit cards sets 6.0 through 13.0 if no equations of motion matrices are to be modified. 

Card sets 6.0 through 13.0 contain operational instructions and data used to modify the 
equations of motion for use in the solution program L221 (TEV156) or any other 
program that is compatible with these output results. 


Card Set 6.0 - Equations of Motion Data 


COLS. 


KEYWORD/ 

VARIABLE 


FORMAT 


DESCRIPTION 


I- iO $EOM 

II- 20 INEOM 


21-25 


INEOMF 


A 4 , 6X 
A7,3X 


15 


Keyword introducing the data for equations of motion 
File name where input equationsof motion matrices reside. 
(Default: INEOM = EOMTAP) 

File position number where equationsof motion matrices 
resides. 


26-65 


DYLOFLEXl 

or 

NULEOM, 


II 


A4 , 36X 
or 
815 


(Default: INEOMF = 1) 

Keyword DYLOFLEX indicates that the null matrix indicator 
array is read from the file on which the equations of motion 


matrices reside. Otherwise, the null matrix indicator array 


is read from this card. 


NULEOMj = 0, matrix is null and omitted from file INEOM. 


NULEOMg / 0, matrix is to be read from file INEOM. 

NULEOMg corresponds to 

NULEOM2 corresponds to M2 

NULEOM ^ corresponds to M3 

NULEOM4 corresponds to M4 

NULEOMg corresponds to Mg 

NULEOMg corresponds to C3 

NULEOM 7 corresponds to fj_ 

NULEOMg corresponds to $ 

(No default) 
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Card Set 7.0 - Stability Derivative Data (Optional) 

Card 7.1 - Derivative (Stability) Data 

If this option is used and the gust zones (NPAN, card set 3.0) = 1 (no gust penetration), 
the forcing function coefficient matrix (NPAN=1) or C3 (NPAN=0) is modified to be 
compatible with the response generalized forces M4 and M5. However, if the gust zones 
> 1 (gust penetration), the forcing function coefficient matrix is not modified and 
should be modified manually by using card sets 11.0 or 12.0 to be consistent with the 
response generalized forces. If the gust coefficient modification is not performed, errors 
may result in the responses and loads of the coordinate. 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

DERIVATIVE 

A4,6X 

Keyword introducing derivative (stability) data. 

11-20 

FROM 

10X 

Descriptive word (preposition) 

21-30 

ICARDl 
ITAPE f 

A4,6X 

Keyword CARD indicates derivatives are to be input on cards 




7.5.1 through 7.7.2. 




Keyword TAPE indicates derivatives are to be input on a 




FLEXSTAB* SDSSTP file (card 7.4.1). 

31-35 

NCS 

15 

Number of control surfaces. If control surface 




derivatives are requested (NCS > 0) read card 7.4.2 for 




control surface names. 




(Default: NCS - 0) Maximum « 20. 

36-40 

INDUN 

15 

Indicator to request unsteady derivatives. 




INDUN = 0, do not read unsteady derivatives 




INDUN / 0, read cards 7. 8. 1-7. 8. 2 for unsteady derivativsa 




(Default: INDUN = 0) 

41-50 

QUEBAR 

E10.0 

Dynamic pressure, q”, (force/length^) ** 




(Default: QUEBAR from EOMTAP ) 

51-60 

VT 

E10.0 

Velocity, true air speed, V T , (length/sec. ) ** 




(Default: VT from EOMTAP) 


♦Indicates the FLEXSTAB program with the DYLOFLEX modifications is incorporated 
into the SD&SS program. 


**The units of force and length must be consistent and identical throughout this 
program and the units of the input matrices. 
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Card 7.2 - Input Matrix Column Numbers of Rigid Body Freedoms 


COLS, 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-5 

flXCOLl* 

\lYCOLl 

15 

Column number of the |yj 

freedoms 




(Default: column element 

not changed) 

6-10 

fIZCOLl* 
flPCOL J 

15 

Column number of the 

freedoms 




(Default: column element 

not changed) 

11-15 

r itcolI* 

1 ISCOL.f 

15 

Column number of the J 

freedoms 




(Default: column element 

not changed) 

16-70 

IDCOLj 

1115 

Column number of the control surface freedom. 



(If more 
cards are 
needed, 
FORMAT 
for cards 
following 
is 1415) 

(1=1, NCS see card 7.1) 
(Default: column element 

not changed) 


^Throughout cards 7.2 through 7.8.1, the upper number in brackets is for symmetric 
analysis, and the lower number is for antisymmetric. 

If any column numbers are left blank, no changes will be made to the matrix elements 
of those freedoms. 



Card 7.3 - Constants Associated with Derivatives 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

AX 

E10.0 

X distance from the stability derivative reference 
point to the rigid body motion reference point, + AFT 
(length)** 

(Default: Ax * 0.0) (See fiqure 4.0) 

11-20 

AZ 

E10.0 

Z distance from the stability derivative reference 
point to rigid body motion reference point, + up 
(length) ** 

(Default: Az * 0.0) (See fiquire 4.0) 

21-30 

ALPHA 1 

E10.0 

1G angle of attack, a, (degrees) 
(Default: ALPHA1 « 0.0) 

31-40 

SW 

E10. 0 

Wing reference area, SW (length 2 )** 
(Defaults : 

If CARD* - Fatal Error 
If TAPE* - Value from SDSSTP) 

41-50 

{CBARj 

E10.0 

Reference chord, "c, (length)** 
Reference span, b, (length)** 
(Defaults : 

If CARD* - Fatal Error 
If TAPE* - Value from SDSSTP) 

51-60 

CL1R 

E10.0 

"RIGID" steady state derivative, C L1R 
(Defaults : 

If CARD* - CLlR * 0 and warning message printed. 
If TAPE* - Value from SDSSTP) 

61-70 

CL1E 

E10.0 

— 

"ELASTIC INCREMENT" steady state derivative, C Ll£ 
(same defaults at C £ ^ R 


*Note: Keyword CARD or TAPE is defined on Card 7.1, cols. 21-30. 
**See note on Card 7.1 
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Ill 


Card 7.4.1 - Files Containing FLEXSTAB Aerodynamic Data 
Read this card only if columns 21 through 30 are TAPE on card 7.1. 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

SDSStp 

A10 

Keyword introducing replacement name for the FLEXSTAB 
aerodynamic data file. 

11-20 

SDINDX 

A10 

Name of the file (tape or disk) containing the index 
table to the FLEXSTAB aerodynamic data. 

The name must begin in column eleven, begin with a 
letter, and have less than seven characters. 

Default: SDINDX = "SDINDX" 

21-30 

SDDATA 

A10 

Name of the file (tape or disk) containing the FLEXSTAB 
aerodynamic data. NOTE: SDINDX / SDDATA. 

The name must begin in column 21, begin with a letter 
and have less then seven characters. 

Default: SDDATA = "SDDATA" 

31-35 

36-70 

ISCAS 

15 

FLEXSTAB aerodynamic data case number from which L217 (EOM) 
must extract derivative. 

Default: ISCAS =1 

Available for comments 


Card 7.4.2 - Derivatives from File per Control Surface 


Read this card if columns 21 through 30 are TAPE, and NCS > 0, on card 7.1. Then go 
to card set 7.8 for unsteady derivatives. 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-70 

NAMESCI 

7A10 

(Repeat 

this 

card if 

more 

names 

are 

needed) 

User defined active control surface names associated 
with active control surface derivatives. 

(1=1, NCS see card 7.1) 

Names are first defined in FLEXSTAB. 
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Note: Read cards 7.5.1 through 7.8.2 only if keyword CARD (card 7.1) is selected. 
Card 7.5.1 - RIGID Derivatives from Cards 


11-20 


21-30 


31-40 


41-50 


51-60 


KEYWORD/ 

VARIABLE 


n 

IcybJ 

I CDU 1 

clbrefJ 

( CMUREFl 
CNBREFj 

( CLA 1 

CYP J 

1 CDA 1 

clprefJ 

( cmarefI 
cnprefJ 


FORMAT 


DESCRIPTION 


E10.0 Steady state derivative 
E10.0 Steady state derivative 
E10 . 0 Steady state derivative 
E10.0 Steady state derivative 
E10.0 Steady state derivative 
E10.0 Steady state derivative 


Card 7.5.2 - ELASTIC INCREMENT Derivatives Read from Cards 
Read this card for SYMMETRIC analysis only (card set 5.0). 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

CLUE 

E10.0 

Steady state derivative, C T ~ 

L E 

11-20 

21-30 

CDUE 

dummy 

E10.0 

Steady state derivative, C d q 

E 

31-40 

CLAE 

E10.0 

Steady state derivative, C^ a 

41-50 

CLAE 

E10.0 

Steady state derivative, C Da 

E 

















Card 7.6.1 - RIGID Derivatives Read from Cards 


COLS. 


KEYWORD/ 

VARIABLE 


FORMAT 


DESCRIPTION 



E10.0 


ElO . 0 


E10.0 


Steady state derivative, 


Steady state derivative. 


Steady state derivative, 


L $R 

y e 

K 

NW 

n REF 


: 


Card 7.6.2 - ELASTIC INCREMENT Derivatives Read from Cards 
Read this card for SYMMETRIC analysis only (card set 5.0). 



KEYWORD/ 



COLS. 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

CLQE 

ElO. 0 

Steady state derivative, C A 

L< *E 

11-20 

CDQE 

E10.0 

Steady state derivative, C 


Card 7.7.1 - RIGID Derivatives Read from Cards per Control Surface 

Repeat card 7.7.1 and 7.7.2 in pairs NCS times (card 7.1.). If NCS = 0, omit this card. 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

(- 

i 

H* 

O 

( CLDi j 

E10.0 


IcYDjj 


11-20 

fCDDi 1 

[cLDREFj J 

E10.0 

21-30 

fCMDREFj! 

IcNDREFj J 

E10.0 


DESCRIPTION 


Steady state derivative. 


Steady state derivative. 


(N 

1 yfi T 


'D6 ] 

: ie 


REF, 


Steady state derivative. 


f-W, 

n<S REF T 


















Card 7.7.2 - ELASTIC INCREMENT Derivatives Read from Cards per Con 
Read this card for SYMMETRIC analysis only (card set 5.0). 


COLS. 

keyword/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

CLDEj 

E10.0 

Steady state derivative, C . 

r t 

11-20 

CDDEj 

E10.0 

Steady state derivative, C , A 

° Rt 




(I * 1, NCS) 1 


Card 7.8.1 - RIGID Unsteady Derivatives Read from Card 
Read this card if INDUN /= 0 (card 7.1). 



Card 7.8.2 - ELASTIC INCREMENT Unsteady Derivatives Read from Car< 
NOTES: 1. Read this card for SYMMETRIC analysis only (card set 5.0) 

2. Read this card if INDUN = 0 (card 7.1). 

KEYWORD/ 

COLS. VARIABLE FORMAT DESCRIPTION 

1-10 CLADTE El 0.0 


11-20 


CDADTE 


E10.0 










Card Set 8.0 - Scale EOM Matrix Elements (Optional) 

This card is repeated for each EOM matrix to be scaled. 


COLS. 

KEYWORD/ 

VARIABLE 

— 

FORMAT 

DESCRIPTION 

1-10 

SCALe 

A10 

Keyword introducing matrix elements to be scaled. 

11-20 

MATNAM 

A10 

Equations of motion matrix name to be scaled. Matrix 




name must be one of the following keywords? 




Ml, M2, M3, FREQ, M4 , M5 , C3 , FL, PHI 

21-25 

IFREQ 

15 

Matrix of Ith frequency 




(Default: IFREQ = 1) 

26-30 

dummy 

5X 

Blanks 

31-40 

SCLMAT 

E10.0 

Scalar by which each element of this matrix is 




multiplied. 


Card Set 9.0 - Sensor Data (Optional) 
Card 9.1 - Introduce Sensor Data 


COLS. 


1-10 

11-20 

21-25 

26-30 

31-45 


KEYWORD/ 

VARIABLE 


SENSor 


INSENF 


NLDSEN 


DYLOFLEX 

or 

NULSEN 


FORMAT 


A4 , 6X 
A7, 3X 

15 

15 


(A4,11X) 

f 315 / 


DESCRIPTION 


Keyword introducing sensor data 

Name of file on which the input sensor matrices reside. 
(Default: INSEN - LODTP2) 

File position number in which sensor matrices reside. 
(Default; INSENF = 1) 

Number of loads (row size) of the sensor matrices being 
read from INSEN file. 

(No default) 

Keyword DYLOFLEX indicates that the null matrix indicator 
array is read from the file on which the sensor matrices 
reside, otherwise, the null matrix indicator array is read 
from this card: 

NULSEN j * 0, matrix is null and omitted from the file. 
NULSENj / 0, matrix is read from file. 

NULSEN ^ corresponds to 

NULSEN ^ corresponds to m 2 

NULSEN ^ corresponds to 
(no default) 
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Card 9.2 - Row Selection of Sensor Data 

Repeat this card for each matrix from which sensor rows are selected. 


COLS . 


KEYWORD/ 

VARIABLE 


FORMAT 


DESCRIPTION 


1-10 


f MlBAr 1 


, M2BAr L 
M3BA r I 


A4,6X 


11-15 


INROW 


15 


• 16-20 


IUTROW 


15 


Matrix name of sensor data. Matrix name must be one of 
the following keywords : 

MlBAR, M2 BAR, M3BAR 
(No default) 

The row number of the input matrix where sensor data will 
be selected. 

(No default) 

The row number of the output (augmented) matrix where the 
selected sensor data will be placed. 

(No default) 


( 21-25 
26-30 

1 31-35 
36-40 

1 41-45 
46-50 

( 51-55 

56-60 

i 

( 61-65 
66-70 


INROW 

IUTROW 

INROW 

IUTROW 

INROW 

IUTROW 

INROW 

IUTROW 

INROW 

IUTROW 


15 

15 

15 

15 

15 

15 

15 

15 

15 

15 


1 


Card 9.3 - End Sensor Data 


Each pair of input and output row numbers is repeated 
for each row selection. 

If more than six rows selected, repeat this card 
with the same matrix name. 

NOTE: a "-1.” is also placed in the "IUTROW" row - 

column diagonal element of the Ml equation of 
motion matrix. 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

END sensor 

A4,6X 

Keyword indicating the end of the sensor data. 
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Card Set 10.0 - SAS Data (Optional) 
Card 10.1 - Introduce SAS Data 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

SAS 

A4,6X 

Keyword introducing SAS data. 


Card 10.2 - SAS Data 

This card is repeated as many times as necessary in order to define the elements of the 
stability augmentation system (SAS) equations which are to be placed in the equations 
of motion. 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-5 

ISAS 

15 

Ith row of augmented equation to the equations of notion 
(Default: Pickup previous row number; exception, the 

first ISAS has no default.) 

6-10 

JSAS 

15 

Jth column of augmented equation to the equationsof motion 
(No default) 

11-20 

Ml IJ 

E 10.0 

Value of the element of the SAS equation to be placed in 
the Ith, Jth location of the M^ matrix. 

21-30 

M2IJ 

E 10.0 

Value of the element of the SAS equation to be placed in 
the Ith, Jth location of the M^ matrix. 

Ul 

1 

Xr- 

O 

M3IJ 

E 10.0 

Value of the element of the SAS equation to be placed in 
the Ith, Jth location of the M 3 matrix. 


Card 10.3 - End SAS Data 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

END sas 

A4 , 6X 

Keyword indicating all SAS data has been defined. 
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Card Set 11.0— Replace EOM Matrix Elements (Optional) 

This card may be repeated. 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

REPLace 

A6,6X 

Keyword introducing matrix elements to be replaced. 

11-20 

MATNAM 

A10 

Eguationsof motion matrix name to be replaced. Matrix 




name must be one of the following keywords: 




Ml, M2, M3, FREQ, M4 , M5, 




C3 , FL, PHI 

21-25 

IFREQ 

15 

Matrix of Ith frequency 




(Default: IFREQ = 1) 

26-30 

dummy 

5X 

Blanks 

31-35 

I ROW 

15 

Ith row of matrix 

36-40 

J COL 

15 

Jth column of matrix 

41-50 

AIJ 

E10.0 

Value replacing Ith, Jth element of matrix 


Card Set 12.0 - Increment EOM Matrix Elements (Optional) 
This card may be repeated. 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

. 1-10 

INCRement 

A4 , 6X 

Keyword introducing matrix elements to be incremented 

11-20 

MATNAM 

A10 

Equation of motion matrix name to be incremented. Matrix 




Name must be one of the following keywords: 




Ml, M2, M3, FREQ, M4 , M5 , 




C3 , FL, PHI 

21-25 

IFREQ 

15 

Matrix of Ith frequency 




(Default: IFREQ = 1) 

26-30 

dummy 

5X 

Blanks 

31-35 

I ROW 

15 

Ith row of matrix 

36-40 

JCOL 

15 

Jth column of matrix 

41-50 

AIJ 

E10.0 

Value incrementing Ith, Jth element of matrix 
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If the user wishes to transform the equations of motion and load equations from inertial 
axes to body-fixed axes, this card set is used. This card set is not repeated. 


Card Set 13.0 - Body Axis Transformation (Optional) 

For this card set, the upper number in brackets is for the symmetric analysis, and the 
lower number is for the antisymmetric analysis. 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

BODYaxis 

A4 ,6X 

Keyword introducing body axis transformation data 

11-15 

flCOLXl 

IicolyJ 

15 

| C0L | column of the original axis 
YcOL J (Default: column elements not used) 

16-20 

flCOLZl 

[ICOLPJ 

15 

I COL J co i umn Q f the original axis 

14>co L J (No defau i fc ) 

21-25 

riCOLTl 

llCOLsJ 

15 

J COL | column of the original axis 

If’coL J (No default > 

26-30 

dummy 

5X 

Blanks 

31-40 

ALPHAl 

E10. 0 

1G angle of attack (degrees) | 




(Default: ALPHAl = 0.0) 

41-50 

BODYVT 

E10.0 



Velocity, true air speed, V T , (length/sec.)* 
(Default: Body V T from EOMTAP ) 


*See note on card 7.1. 
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6.3.3 INSTRUCTIONS TO MODIFY LOADS MATRICES 

Omit card sets 14.0 through 17.0 if no load equation matrices are to be modified. 

Card sets 14.0 through 17.0 contain operational instructions and data used to modify 
the load equations for use in the solution program L221 (TEV156) or any other program 
that is compatible with these output results. 

Card Set 14.0 - LOADS Equations Data 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

1-10 

$LOAds 

A4 , 6X 

11-20 

INLOD 

A 7 , 3X 

21-25 

INLODF 

15 

26-30 

NLDOU 

15 

31 ■'•65 

f DYLOFLX 1 
iNULLODj f 

IA4,3X\ 
1715 / 


DESCRIPTION 


Keyword introducing the data for load equations 
File name where input load equationsmatrices reside. 
(Default: INLOD = LODTAP) 

File position number where load equationsmatrices reside 
(Default: INLODF = I) 

Total number of output loads 
(No default) 

Keyword DYLOFLX indicates that the null matrix indicator 


array is 

read from the file 

on which the load 

equations 

jnatrices 

reside. Otherwise, 

, the null matrix 

indicator 

array is 

read from the 

same 

card columns. 


NULLODj 

= 0, matrix is 

nul 1 

and omitted from 

the file. 

NULLODj 

/ 0, matrix is 

read 

from file. 


NULLOD.J 

corresponds to 

“l 



nullod 2 

corresponds to 

r 2 



nullod 3 

corresponds to 

*3 



NULLOD y 

corresponds to 




nullod 5 

corresponds to 

^5 



NULLODg 

corresponds to 

C 3 



NULLOD ? 

corresponds to 

* 




Note: Card set 14.0 is required if SAS equations are added to the equations of motion. 
The use of card set 14.0 will increase the column size of the load coefficient 
matrices. 
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Card Set 15.0 - SCALE LOADS Matrix Elements (Optional) 


Repeat this card for each LOADS matrix to be scaled. 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

SCALe 

A10 

Keyword introducing matrix elements to be scaled. 

11-20 

MATNAM 

Al 0 

Load equation matrix name to be scaled, matrix name must 




be one of the following keywords: 




Ml BAR, M2BAR, M3BAR, M4BAR , M5BAR , C3BAR , PHIBAR 

20-25 

IFREQ 

15 

Matrix of Ith frequency 




(Default: IFREQ = 1) 

26-30 

dummy 

5X 

Blanks 

31-40 

SCLMAT 

El 0. 0 

Scalar, multiply each element of this matrix by this value. 


Card Set 16.0 - Replace LOADS Matrix Elements (Optional) 
This card may be repeated. 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

REPLace 

A4 , 6X 

Keyword introducing matrix elements to be replaced. 

11-20 

MATNAM 

Al 0 

Loads Equation matrix name to be replaced. Matrix name 




must be one of the following keywords: 




M1BAR, M2 BAR , M3BAR , M4BAR , M5BAR, C3BAR , PHIBAR 

21-25 

IFREQ 

15 

Matrix of Ith frequency 




(Default: IFREQ = 1) 

26-30 

dummy 

5X 

Blanks 

31-35 

I ROW 

15 

Ith row of matrix 

36-40 

JCOL 

15 

Jth column of matrix 

41-50 

AIJ 

E10.0 

Value replacing Ith, Jth element of matrix 
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Card Set 17.0 - Increment LOADS Matrix Elements (Optional) 

This card may be repeated. 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

INCRement 

A4 , 6X 

Keyword introducing matrix elements to be incremented. 

11-20 

MATNAM 

A10 

Load equation matrix name to be incremented matrix name 




must be one of the following keywords : 




M1BAR, M2BAR, M3BAR , M4BAR , M5BAR, C3BAR , PHIBAR 

.21-25 

IFREQ 

15 

Matrix of Ith frequency 




(Default: IFREQ = 1) 

26-30 

dummy 

5X 

Blanks 

31-35 

I ROW 

15 

Ith row of matrix 

36-40 

JCOL 

15 

Jth column of matrix 

7 

s 

AIJ 

£10.0 

Value incrementing Ith, Jth element of matrix 


6.3.4 INSTRUCTIONS FOR PREPARATION OF QR MATRICES 

Omit card set 18.0 if no QR matrices are to be generated. 

Card set 18.0 contains operational instructions and data used to modify the equations of 
motion and load equations for use in the solution program QR or any other program 
that is compatible with these output results. 

Card Set 18.0 - QR Data Preparation (Optional) 

Repeat cards 18.1 and 18.2 for each of either cards 18.3 or 18.4 or 18.5, and 18.6 

1. QR data may be executed with the modified or unmodified EOM and LOADS 
matrices as specified on cards 18.1 and 18.6 respectively. 

2. QR data may be executed without reading card sets 4.0 through 17.0 
Card 18.1 - Request QR Data Preparation 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

$QR 

A4 , 6X 

Keyword introducing the data for QR output 

11-20 

IUTQR 

A7, 3X 

File name where QR matrices are to be written. 




(Default: IUTQR = QRTAP) 

21-25 

IFLQR 

15 

File position number where QR matrices are to be 




written . 




(Default: IFLQR = 1) 


53 



Card 18.2 - Input Equations of Motion Matrices for QR 


COLS. 


1-10 


11-20 


21-25 


26-65 


KEYWORD/ 

VARIBLE 


QREO m 


INEOM 


INEOMF 


/dyloflx! 

INULEOMj I 


FORMAT 


A4 , 6X 


A7 , 3X 


15 


/A4,36xl 
\8I5 f 


DESCRIPTION 


Keyword introducing the source of the equationsof 
motion matrices. 

File name where input equationsof motion matrices reside. 
{Default: INEOM = EQEOM) 

File position number where equationsof motion matrices 
resides 

(Default: INEOMF =1) 

Keyword DYLOFLX indicate null matrix indicator array is read 
from file where equationsof motion matrices reside. Other- 
wise, null matrix indication array is read from cards. 


NULEOM, 


= 0, matrix is null and omitted from the file. 


NULEOM^ / 0, matrix is read from file 

NULEOM^ corresponds to 

NULEOM^ corresponds to M 2 

NULEOM ^ corresponds to M 3 

NULEOM^ corresponds to 

NULEOM,. corresponds to M,. 

NULEOM^ corresponds to 
NULEOM^ corresponds to f^ 

NULEOM corresponds to f 

O 

(No default) 
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Card 18.3 - QR Wagner Option 


This card is used if roots with Wagner indical lift growth functions are to be calculated 
by QR. 


COLS . 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

; i-io 

WAGN er 

A4 , 6X 

Keyword indicating QR matrices are to be formulated using the 




equations of motion with Wagner functions. 

11-20 

AQR 

E10.0 

Wagner function, aj 

21-30 

BQR 

E10.0 

Wagner function, b^ 




(see equation 15 ) 

31-40 

ALQR 

E10.0 

Wagner function, 

;41-50 

BEQR 

E10.0 

Wagner function, 8^ 


Card 18.4 - QR Root Option 

This card is used if roots without Wagner indical lift growth function are to be 
calculated by QR. 


COLS. 

keyword/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

ROOT 

A4,6X 

Keyword indicating 

QR matrices are to be formulated using the 




equations of motion 

without Wagner functions. 


Card 18.5 - QR TIME Option 

This card is used if time histories are to be calculated by QR. 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

TIME 

A4,6X 

Keyword indicating QR matrices are to be formulated using the 
equations of motion and load equation matrices for a time 
history solution. 
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Card 18.6 - Input Load Equations Matrices for QR 
This card is required if keyword TIME (card 18.5) is used. 


COLS. 


31-65 


KEYWORD/ 

VARABLE 


1-10 


11-20 


21-25 


26-30 


QRLOad 


INLOD 


INLODF 


NLDQR 


DYLOFLX 

NULLOD, 


FORMAT 


A4 , 6X 


A7 , 3X 


15 


15 


Att , 3X 
715 


DESCRIPTION 


Keyword introducing the source of the load 
equation matrices. 

File name where input load equation matrices reside 
(Default: INLOD - EQLOD) 

File position number where load equationsmatrices reside 
(Default: INLODF - 1) 

Number of loads 

(Default: If DYLOFLX file, NLDQR is extracted from the first 

recordon the file INLOD). 

Keyword DYLOFLX indicate null matrix indicator array is 
read from file - where load equationsmatrices reside. 

Otherwise, null matrix indicator array is read from this card. 


NULLODj 

- 0, matrix is 

null and omitted from the file 

NULLODj 

W 0, matrix is 

to be read from file INI/3D 

NULLOD^ 

corresponds 

to 

M 1 

nullod 2 

corresponds 

to 

M 2 

nullod 3 

corresponds 

to 

M 3 

NULLOD^ 

corresponds 

to 

M « 

nullod 5 

corresponds 

to 

M 5 

NULLODg 

corresponds 

to 

C 3 

NULLOD ? 

corresponds 

to 

* 


Card Set 19.0 - Terminator 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

$QUI t 

AM , 6X 

Keyword indicating the last data of the EQMOD module has 
been read. The $QUIT card may include comments following 
the fourth column. 
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6.3.5 SUMMARY OF CARD INPUT DATA 
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Requi rements Reference 

or Function Key Words and/or Data Card Format Card Set (CS) 
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6.4 MAGNETIC FILES INPUT DATA 


The input matrices to the L219 (EQMOD) program will normally be obtained from 
magnetic files (tape or disk) prepared by the programs; Equations of Motion, L217 
(EOM), and Load Equations, L218 (LOADS). A magnetic file prepared by FLEXSTAB 
with the DYLOFLEX modifications will also be needed when using the stability 
derivatives generated in FLEXSTAB. However, because the EQMOD output magnetic 
files are in the same format as the input files, it is possible for EQMOD to use as input 
the magnetic files generated by a previous execution of EQMOD. In addition, any user 
generated magnetic files(s) may be used as input into EQMOD if they have the required 
format. 

The format for the equations of motion input magnetic file is shown in figure 8. The 
format for the load equations input magnetic files is shown in figure 9 for the load 
equations, and in figure 10 for the sensor equations (which are themselves a specific 
type of load equation). 

All EOM and LOADS matrices are in the READTP/WRTETP format. 1 

The stability derivatives generated in FLEXSTAB are written on the SDSSTP magnetic 
file. Before executing L219 (EQMOD), this data must be copied by EQMOD onto two 
magnetic files, SDINDX and SDDATA. SDINDX contains the FLEXSTAB index matrix 
and SDDATA the FLEXSTAB stability derivative data shown in figure 11. 

6.5 OUTPUT DATA 


6.5.1 PRINTED 

All card input data will be printed as read and interpreted. Optionally the following 
data may be printed: 

1. Matrices read from files 

2. Matrices written on output files. 

6.5.2 MAGNETIC FILES 

EQMOD will write as many as three magnetic files containing modified equations of 
motion, modified load equations, and a file for the Linear System Analysis program, 
QR, which may contain the characteristic equations for the modified equations of 
motion and/or the modified equations of motion and load equations. 


^R. E. Clemmons: Programming Specifications for Modules of the Dynamic Loads Analysis System 
to Interface with FLEXSTAB . NASA contract NAS1-13918, BCS-G0701 (internal document) , 
September 1975. 
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SIZE 



NFREQH - Number of frequencies 
NPANI - Number of panels input 

* Note; If {freqm} exists then C 5 and will be complex 


Figure 8. - Equations of Motion Input Fite ( EOMTAP ) 
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SIZE 

30 x 1 

NLDI x NDOFI 
NLDI x NDOFI 
NLDI x NDOFI 
NLDI x NDOFI 
NLDI x NDOFI 

NLDI x 1 
NLDI x NPANI 


where; 



NIDI - Number of loads input 

NDOFI - Number of degrees of freedom input 

NPANI - Number of panels input 

NFREQM - Number of frequencies 

NLOD - Number of load sets 


* Note: If NFREQM > 1 then C* and $ will be complex 
Figure 9. — Load Equations input Fife (LODTAP) 


Repeated 
for each 
load set 
(NLOD) 
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30 x 1 


NSENS 

NSENS 

NSENS 


where : 


x NDOFI 


{Header} 

[Mi] 


x NDOFI 
x NDOFI 



Repeated 
for each 
load set 

(nlod) 


NSENS - Number of sensor loads 

NDOFI - Number of degrees of freedom input 

NLOD - Number of load sets 


Figure 10. — Sensor Equations Input File (L0DTP2) from the A VD Loads Path 


FLEXSTAB 

name 


Engineering symbol 


fCONTRL)— A 


(CONTRL)-S 


ICSNAMES)A 


(CSNAMES)S 


c v. 

"rigid 

S, . 

elastic 

c \ 

rigid 

Cy 

6 r 

elastic 


incr 


incr 

c e 

V ■ 

a rigid 

\ 

elastic 

Cn 

rigid 

Cp 

elastic 


incr 


incr 

\ . . 
a rigid 

\ . 
elastic 

C 

5 r 

rigid 

S> 

elastic 


incr 


incr 

where i ranges over the number of antisymmetric active contr 

V. 

rigid 

c i 

L « e . 
elastic 

incr 

C l 

% 

rigid 

Ci 

% 

elastic 

incr 

c n 

5 e. . 
e rigid 

X . 

elastic 

incr 

c \ 

'rigid 

C D 

s 

elastic 



incr 

C m 

"V . 

rigid 

C m 5 

e etastic 

incr 

c 

mt 

5 Cj 

rigid 

s 

elastic 



incr 


*0, . . X 

rigid elastic 

incr 


c c- 

'rigid 'elastic 


C_ C 

\ n «c 

rigid 


incr 

'elastic 

incr 


where i ranges over the number of symmetric active controls, i=l 
JTlLERON RUDDER AC. 


'namej name 


n 


where i ranges over the number of antisymmetric active controls, i-1 ,n _ 2 


El 


ELEVATOR AC n . m AC 

namej name 


a 


where i ranges over the number of symmetric active controls, i 52 ! ,n - 1 


Figure 1 1. — Genera i Form of Derivative Matrices on Input File (SDSSTP) 
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FLEXSTAB 






name 


Engineering symbol 



(STATIC)-S 


F l - ■ •„ 

rigid 

Ci Ci A Ct A 

*-i L u L i i 

elastic rigid elastic 

C L 

0 rigid 

C L c lA C lA 

°elastic ^ rigid ^elastic 





incr incr 


incr 

incr 




0 

0 C D j C Df) 

u rigid u elastic 

C D q 

rigid 

C D C DA C D^ 

“elastic p rigid ^elastic 





incr 


incr. 

incr 




0 

0 C A C A 

mn m y 

rigid elastic 

c m 

m 

“rigid 

C C A C A 

m Q mJL' 

“elastic p rigid ^elastic 





incr 


incr. 

incr 








— 


{STATIC)-A 


frigid 

C c yA 

^elastic p rig»d 

CyA 

p e1astic 

Cy A 

r rigid 

C yf v 

elastic 





incr 

incr 


incr 




Co 

frigid 

C^ C fA 

^elastic p rigid 

C gA 

p elastic 

^j?A 

rigid 

C 2A 

r elastic ' 





incr 

incr 


incr 




Sigid 

C n e C nA 

^elastic ^ rigid 

C nA 

p elastic 

C nA 

rigid 

elastic 




— 

incr 

incr 


incr 



Figure 1 1. — (Concluded) 





The files of the modified equations of motion and loads equations have the same format 
as the input formats shown in figures 8 and 9 except that the size of the modified 
matrices may be larger by: 

NDOF =NDOFI + NDOF NSEN + NDOF NS AS 


where: 


NDOF 

NDOFI 

NDOFnseN 

ndofnsas 


= The total number of degrees of freedom 
= The number of degrees of freedom input 
= The sensors number of degrees of freedom 
— The SAS number degrees of freedom 


The default name is EQEOM for the modified equations of motion file and EQLOD for 
the modified load equations file. 


'I'he QR equations output magnetic file format is shown in figure 12. 

AU matrices contained on these output magnetic files are written by the WRTETP 
subroutine. 


6.6 RESTRICTIONS 

The following restrictions apply only when generating matrices for QR. These 
restrictions are due to the limitation of the QR program (ref. 6). 

• No gradual penetration is allowed when forming matrices for the time history 
solution in QR. 

No multi-forcing-function is allowed for a QR time history solution. 

• Only first frequency used (EQMOD restriction). 

• Only one load set used (specified via card). 

6.7 DIAGNOSTICS 

6.7.1 FATAL ERRORS 

All fatal errors detected by L219 (EQMOD) will result in the printing of a diagnostic 
error message. These, messages are self-explanatory and are of the following format: 

******** FATAL ERROR (nnnnn) DIAGNOSED WHILE EXECUTING ROUTINE 
(routine name). 

Any additional error message follows. 
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SIZE 


NDOF x NDOF 

[S-] 


NDOF x NDOF 

[S 3 ] 


NDOF x NDOF 

[s 2 ! 


NDOF x NDOF 

[s 5 ] 


NDOF x NDOF 

[s°] 



End-of-Fi 1 e 


NDOF x NDOF 

[s 2 ] 


NDOF x NDOF 

[S'] 

> 

NDOF x NDOF 

[s :] 



End-of-Fi 1 e 




' 

NDOFL x NDOFL 

[s 2 ] 


NDOFL x NDOFL 

[S'] 

> 

NDOFL x NDOFL 

[s°] 


NDOFL x 1 

[C>] 



End -of- Fi 1 e 


— 




File generated 
for Wagner option 


File generated 
for Root option 


File generated 
for Time option 


where : 

NDOF - Number of degrees of freedom 

NDOFL - NDOF + NLD 

NLD - Number of loads 


Figure 12. — QR Equations Output File (QRTAP) 
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where nimnn is a diagnostic number from the following list: 


Error 

number 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 


Description 

$EQMOD is not the first card input data 

Error returned from FETAD (see FETAD error message) 

Keyword (aaaaaaaaaa) with code number (nnnnn) is not recognized. 

Total number of degrees of freedom (nnnnn) not within range of 1 and 100. 
Total number of panels (nnnnn) not within range of 0 and 50. 

Number of frequencies (nnnnn) not within range of 0 and 20. 

Output file name (aaaaaaaaaa) for modified EOM matrices, invalid. 

Output file name (aaaaaaaaaa) for modified LOADS matrices, invalid. 

The requested number of (LOADS/EOM/QR) sets (nnnnn + 1) is greater 
than the maximum allowed (nnnnn). 

READTP error number (nnnnn) returned. (See section 6.7.3.) 

Dimension on first matrix of DYLOFLEX (EOM/LOADS) input tape is not 
30 x 1. 

(EOM/LOADS) scale data must be grouped together. 

(EOM/LOADS) replacement data must be grouped together. 

(EOM/LOADS) increment data must be grouped together. 

Body axis data already defined. 

Keyword (aaaaaaaaaa) is not recognized as a SENSOR data. 

Sensor matrix (aaaaa) not grouped in increasing (M1BAR,M2BAR, M3 BAR) 
order. 

(keyword) matrix name (aaaaa) not recognized; (keyword) matrix name 
(aaaaa) and frequency number (nnnnn) must be grouped together. 

QR option not specified; must have one of the following (WAGNER, 
ROOT, TIME). 

Field length requested is too small for this problem. 

WRTETP error number (nnnnn) returned. (See section 6.7.4.) 

Matrix (aaaaaaaaaa) frequency (nnnnn) is greater than number of 
frequencies (nnnnn). 

Derivatives input; VT, B, or CBAR contain zero. 

Error returned by FETCHM while reading matrix (aaaaaaaaaa). 

Control surface name specified (aaaaaaaaaa) does not equal to names read 
from SDSSTP (list of names follow). 

QR time option; PHI-TILDA OR PHI-TILDA-BAR (aaaaaaaaaa) has more 
than one column. 
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6.7.2 WARNING MESSAGES 

All warning messages will be self- explanatory and printed in the following format: 


******** WARNING MESSAGE (nnnnn) DIAGNOSED WHILE EXECUTING 
ROUTINE (routine name). 

(Any additional warning message follows.) 

where nnnnn is a warning message number from the following list. 

Warning 

number Description 


1 The maximum number of title cards (nnnnn) has already been read. The 
above title card is treated as a comment. 

2 End of record encountered before $QUIT card. $QUIT card assumed. 

6.7.3 READTP ERROR CODES 


Code 


Description 


0 

2 

3 

4 

5 

6 
7 

1000+1 

3000+1 


No errors detected. 

File or matrix skipping code is negative. 

Dimensioned row size s 0. 

Dimensioned row size is less than the actual row size. 

Name specified does not match the matrix name on the file. 
NROWS*NCOLS ^ 0 or greater than the buffer limit of 10000. 

End-of-file read instead of the two-record matrix. 

Error discovered trying to skip files. I files remained to be skipped when an 
end-of- information (EOI) was encountered. 

Error discovered trying to skip matrices. 1/2 matrices remained to be 
skipped when an end-of-file (EOF) was encountered. 


6.7.4 WRTETP ERROR CODES 

Code Description 


0 

2 

3 

4 
6 

1000+1 

3000+1 


No error detected. 

File on matrix skipping code is negative. 

Dimensioned row size 0. 

Dimensioned row size is less than the actual row size. 

NROWS*NCOLS « 0. 

End-of-information (EOI) was read trying to skip files. I files remained to 
be skipped. 

End-of-file (EOF) was read trying to skip matrices. 1/2 matrices remained to 
be skipped. 
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7.0 SAMPLE PROBLEM 


The sample problem in this section consists of three small test cases (numbers 2, 3 and 
5) formed to exercise most of the options available in the program. The size of these test 
cases allowed them to be easily checked by hand calculations. 

Test Case 2: 

This sample problem assumes that the matrix coefficients for the equations of motion 
were generated by L217 (EOM) for a six degree of freedom system. The load equation 
coefficients are assumed to be read from a file LODTAP that was not generated by the 
DYLOFLEX system. Modifications to the equations include overwriting the rigid body 
terms in the equations of motion with stability derivative data, adding sensor equations 
to the system, and replacing and incrementing matrix elements. Wagner indical lift 
growth effects shall be included, and a data tape for QR rooting and time history 
analysis shall be generated. 

Test Case 3i 

This sample problem assumes ail load and equations of motion coefficient matrices were 
generated by DYLOFLEX programs for a four degree of freedom system. Modifications 
include a scaler multiplication of M 4 , a transformation of the equations into the 
body-fixed axes, and replacement and incrementation of certain load matrix coefficients. 

Test Case 5: 

The last sample problem deals with a seven degree of freedom system. Here the 
equations of motion matrix coefficients will be modified by overwriting the rigid body 
elements with stability derivative data. Sensor and stability augmentation equations 
will also be added. 

. f 


Boeing Commercial Airplane Company 
P.O. Box 3707 

Seattle, Washington 98124 
May 1977 
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SEQH00 

TEST CASE 2 






WO 

TITLE 

CHECK 

CASE 2 , MODIFY DERIVATIVES FROM 

CAROS, ADO 

SENSORS, QR 

TP 

C 

6-OFGR EES OF FREEDOM , 2-FREQUE NCIES 




2 . 1A 

SIZE 

6 

0 2 






3.0 

OUTPUT 

EOEOM 2 

1 

E0L0D2 

l 




4 . 1 

PR INT 

INPUT 

MATRIX 

FREQUENCY 2 




4.2 

PR INT 

OUTPUT 

MA TP ICES 

CHANGED 





4.3 

SYMMETRIC 








5.0 

C 

EQ. OF 

MOTION ASSUMED TO HAVE BEEN GENERATED BY 

L217I EOM) 

2 . IB 

$E0« 

E OMTAP 

2DYL0FLX 





6.0 

C 

ONE CONTROL SURFACE 

AND LSE 

UNSTEADY 

DERI VATI VES 



2. 1C 

DERIVATIVES FROM 

CARD 

1 

i 2.0 

10.0 



7.1 


1 2 

4 






7.2 



. L 

2. 5 

1.0 

3. 0 



7.3 

4.0 

1.0 

1.0 

7.0 

1.0 

10. 0 



7. 5. IP 

1.0 

. 1 

-.1 

-WO 

WO 

-2.0 



7.5.2E 

z.o 

3.0 

7.0 






7 .6.1R 

. 1 

-1.0 

1.0 






7.6.2E 

4.0 

WO 

6-0 






7.7.1R 

-2.0 

- 1.0 

-WO 






7.7.2E 

2.0 

3.0 

5.0 






7.8. IR 

0.0 

0.0 

0.0 






7.0.2E 

C 

SENSOR 

DATA FROM L2 1 8( LOADS) 






SENSOR 

LOOTP 2 

2 

3DYL0FLX 





9.1 

M3BAR 

1 

5 2 

6 





9.2 

END SENSOR 







9.3 

REPLACE 

Ml 


3 

3 200. 




i WO 

INCREMENT 

M3 


3 

3 10. 




12.0-1 

INCREMENT 

M3 


3 

4 5. 




12.0-2 

C 

load tape generated 

OUTSIDE 

OF DYLOFLEX SYSTEM 



2 . IE 

sloads 

lodtap 

1 

2 0 

0 3 

4 5 

8 

0 

14.0 

SCALE 

M4BAR 

2 

2.0 





15.0 

SOR 

OR TAP 

l 






18. IA 

QREOH 

EQE0H2 

DYIOFLX 





18. 2A 

ROOT 








18.4 

SQR 

QRTAP 

2 






18. IB 

OREOM 

E0E0M2 

dyloflx 





18. 2B 

T IHE 








18.5 

QRLOAO 

E0L002 


OYLOFLX 





18.6 

SQR 

QRTAP 

3 






Ifi.lC 

QREOH 

EQE0H2 

dyloflx 





18. 2C 

WAGNER 

.5 

.5 

1.0 

2.0 




18.3 


SQUIT 19.0 
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SEQHOD 

TEST CASE 

3 



TITLE 

CHECK CASE 

3. 

MODIFY 

EOM AND 

SIZE 

4 


2 


OUTPUT 

TESTA 


l 

DUMMY 

PAINT 

OUTPUT 

MATRICES 

ALL 

SYMMETRIC 





SEOM 

EOMTAP 


2DYLOFLX 

SCALE 

M4 



2.0 

BODY AXIS 

1 2 


3 

1.0 

SLOAOS 

LODTAP 


2 

2DYL0FLX 

REPLACE 

M3BAR 



1 

REPLACE 

M4BAR 


2 

1 

INCREMENT 

(QUIT 

M3BAR 



2 



1.0 

LOADS 

3.0 

1 

4.1 


4.3 


5.0 


6.0 


8.0 

10.0 

13.0 

14.0 

4 10.0 

L6.0-1 

4 5.0 

16.0-2 

4 3.0 

17.0 

19.0 



tEQMOO 

TEST CASE 5 




1 .0 

title 

CHECK 

CASE 5, MODIFY DERIVATIVES FROM 

CARDS ANO ADD SAS 

AND SENSORS 

SIZE 

7 

2 




3.0 

OUTPUT 

EQEGM 

1 

EOLOD 

l 


4-1 

PRINT 

INPUT 

MATRIX 

ALL 



4.2 

PR INT 

OUTPUT 

MATR IX 

ALL 



4.3 

SYMMETRIC 





5.0 

$EQM 

EOMTAP 

1 

1 0 

1 1 

110 0 

6.0 

DERIVATIVE FROM 

cards 

1 

0 2.0 

10.0 

7.1 


1 2 

4 




7.2 



• 1 

2.5 

1.0 

3.0 .2 

7.3 

4.0 

1.0 

1.0 

7.0 

1.0 

10.0 

7. 5. IP 

L.O 

. 1 

-.1 

-1.0 

1.0 

-2.0 

7.5.2E 

2.0 

3.0 

7.0 




7.6. 1R 

.1 

-1.0 

1.0 




7.6.2E 

4.0 

1.0 

6.0 




7.7.1P 

-2.0 

-1.0 

-1.0 




7.7. 2E 

SENSOR 



3 0 

1 1 


9.1 

M2BAR 

3 

5 




9.2-1 

M30AR 

1 

6 




9.2-2 

ENO SENSOR 





9.3 

SAS 






10. 1 

7 

7 1-0 





10.2-1 

7 

5 

-1.0 




10.2-2 

7 

6 


2.0 



10.2-3 

4 

7 1.0 

2.0 

3.0 



10.2-4 


ENO SAS 10*3 

$QUI T 


A 




**************************************************** 


PROGRAM L219A2 VERSION JUNE 29,77 NOW RUNNING. 
THE PROGRAM IS PART OF THE OYLOFLX SYSTEM 
DEVELOPED FOR NASA UNDER CONTRACT NASI-13918. 
DATE OF RUN IS 77/11/09. 

TIME OF RUN IS 09.46.16. 


>************«*************************************** 


( SEOMOD 

(TITLE 

(C 

(SIZE 


TEST CASE 2 

CHECK CASE 2, MODIFY DERIVATIVES FROM 
6-DEGREES OF FREEDOM t 2-FREQUENC I ES 
6 0 2 


1.0 

CAROS* ADO SENSORS* QR TP 

2. IA 
3.0 


I 

I 

I 

I 


PROBLEM SIZE 

NDOF * 6* TOTAL NUMBER OF DEGREES OF FREEDOM. 

NP AN - 0* TOTAL NUMBER OF PANELS. 

NFREQM * 2* NUMBER OF FREQUENCIES. 

(OUTPUT EQE0M2 1 EQL0D2 1 


4.1 I 


(PRIM 


OUTPUT TAPE 
IUTEOM * EQE0M2 • 

IFLEOM - l, 

IUTLOO * EQL0D2 * 

IFLLOD * 1* 

INPUT MATRIX 


S 

TAPE NAME FOR MODIFIED EOM MATRICES. 
FILE POSITION NUMBER OF IUTEOM. 

TAPE NAME FOR MODIFIEO LOADS MATRICES. 
FILE POSITION NUMBER OF IUTLOD. 
FREQUENCY 2 


4.2 > 


INPUT MATRICES PRINT OPTION. 

INPR * 2, IF INPR * -999. PRINT ALL INPUT MATRICES. 

* 0, NO INPUT MATRICES PRINTED. 

=* N. MATRICES OF NTH FREQUENCY ONLY PRINTED. 

(PRINT OUTPUT MATRICES CHANGED 4.3 ) 


OUTPUT MATRICES PRINT OPTION. 

IUTPR * 999, IF IUTPR * -999, PRINT ALL OUTPUT MATRICES. 

* 0* NO OUTPUT MATRICES PRINTED. 

* N, MATRICES OF NTH FREQUENCY ONLY PRINTED. 

* 999, PRINT ONLY MODIFIED MATRICES. 


( SYMMETR 1 C 

(C EQ. OF MOTION ASSUMED TO HAVE BEEN GENERATED BY L217IE0M) 

IAEOM EOMTAP 2DYL0FLX 


5.0 
2. IB 

6.0 


I 

) 

I 


EQUATIONS OF MOTION 


IN EOM * EOMTAP , EOM INPUT TAPE NAME. 

INEOMF » 2* FILE POSITION NUMBER OF INEOM. 


(C ONE CONTROL SURFACE AND USE UNSTEAOY DERIVATIVES 2. 1C » 

(DERIVATIVES FROM CARD 1 l 2.0 10.0 7.1 1 


DERIVATIVES FOR SYMMETRIC ANALYSIS. 

I VOL » CARD, INPUT VOLUME 

NCS * 1* NUMBER OF CONTROL SURFACE 

INOUN * 1, UNSTEADY DERIVATIVE INDICATOR 

QUEBAR « . 200E+0 1, DYNAMIC PRESSURE 

VT - • 1 00E+O2, VELOCITY (TRUE AIR SPEED! 


COLUMN NUMBERS OF RIGID BOOY FREEDOMS 
IXCOL * 0, COLUMN OF X FREEOQMS 

1 Z COL - 1* COLUMN OF Z FREEDOMS 

ITCOL * 2, COLUMN OF THETA FREEDOMS 

COLUMN OF DELTA CONTROL SURFACE FREEDOMS FOLLOW 


4 
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CONSTANTS ASSOCIATED WITH DERIVATIVES 

XREF = 0. , X-COORD I NATE OF MOMENT REFERENCE 

ZREF * C. , Z-COORDINATE OF MOMENT REFERENCE 

ALPHA! = • L 75E-02, IG ANGLE OF ATTACK 

SW » .250E+C1, WING REFERENCE AREA 

CBAR » • IQ0E+01, REFERENCE CHORD 

CL 1R = .3C0E+01, RIGID STEADY STATE DERIVATIVE 

CL IE » 0. t ELASTIC STEADY STATE DERIVATIVE 

STEADY STATE DERIVATIVES FROM CARD 
CLU - • 4006 + 01* C— L— LiHA T— P I G ( D 

CDU =» » 1 006 + 0 1 f C-D-UHAT-RIGIO 

CMUREF • • l OOE + O 1 t C- M- UHA T-RI GI O-REF 

CLA = .700E+01, C-L-ALPHA-PIGID 

CDA * ■ 1 OOE+Ol , C-D-ALPHA-RIGID 

CMAREF * . 100E+O2, C-M-ALPHA-RIGID-REF 

CLUE « .ICOE+Ol, C-L— -UHA T- ELASTIC 

CDUE * - l OOE+OO, C-D-UHAT-ELASTIC 

CLAE * -• lOOE+O Lf C-L-ALPHA-ELA STI C 

COAE = . lOOE+Ol, C-D-ALPHA-ELASTIC 

CLQ = . 200E+0 1, C-L-QU6HAT-R1GID 

CDQ » . 3Q0E+OI, C-D-QUEHAT-RIGID 

CMOREF = .700E+01, C-M-QUEHAT-RI GID-REF 

CLOE * -lOOE+OO, C-L-QHAT-ELASTIC 

CDQE = -.100E + 01, C-D— QHAT-ELASTI C 

CONTROL SURFACE DERIVATIVES 
CLO I C-L-DELTA-R IGIDJ FOLLOWS 
.AOOE+Ol 

CDO ( C-O-DELTA-R IGID) FOLLOWS 
. lOOE+Ol 

CMDREF I C-M-DEL TA-R6F I FOLLOWS 
.600E+01 

CLOE IC-L- DELTA— ELASTIC) FOLLOWS 
-.20CE+01 

CDOE I C-C-DELTA-ELA9TIC J FOLLOWS 
lOOE + Ol 

UNSTEADY OER IVAT I VES 

CLADOT = .2000+01, C-L— ALP HA -DOT-HA T-RI G I D 

CDADOT = • 3 QOE+O 1 , C-D- ALPHA -DOT— HA T-RI G I D 

CM ADR F = .5006+01, C- M- ALPHA-DO T-HA T-R I G I D-REF 

CLADTE = 0. , C-L-ALPHA-DOT-HAT-ELASTIC 

COADTE = 0. , C-D-ALPHA-DOT-HAT-ELASTIC 

IC SENSOR DATA FROM L218ILOADSI 

(SENSOR L0DTP2 2 3DYL0FLX 


INSEN 
INSEN F 
NLDSEN 

MATRIX 
(M3BAR 
I END SENSOR 
(REPLACE MI 
( INCREMENT M3 
( INCREMENT M3 


IN OUT 


, SENSORS INPUT TAPE NAME. 

2, FILE POSITION NUMBER OF INSEN. 

3, NUMBER OF LOADS ON SENSOR TAPE. 


3 2 00. 
3 10. 

A 5. 


(C 

( $LQADS 


LOAD tape generated outside of dyloflex system 
LODTAP 1 2 0 0 3 A 5 


LOADS EQUATION 
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INLC1D = LOOTAP 
INLODF = 

NLDOU = 


o LOAOS INPUT TAPE NAME. 

Is FILE POSITION NUMBER CF 1NLCD. 
2s NUMBER OF OUTPUT LOADS- 


(SCALE M4BAR 2 2-0 

( SQR QRTAP t 


15.0 | 

19-14 ) 


0 R 


tUTQR = QRTAP 
IFLQR = 


s TAPE NAME FOR GR MATRICES. 

Is FILE POSITION NUMBER QF IUTQR. 


( QR E CM 

E0E0M2 

OYLOFLX 

1 ROOT 




18. 2A ) 
IB. 4 ) 


OR ROOT OPTION. 
(SOR QRTAP 2 


18. IB ) 


Q R 


IUTQR « QRTAP s TAPE NAME FOR QP MATRICES. 
IFLQR a 2 o FILE POSITION NUMBER OF IUTQR. 


(ORE CM 

E0E0M2 

OYLOFLX 

(TIME 




18.20 ) 
18.5 ) 


Q ft TIME OPTION. 
IQRLQAD EQL002 

(*QR ORTAP 3 


OYLOFLX 


18.6 1 
18. 1C J 


0 R 


IUTQR a QRTAP . TAPE NAME FOR QR MATRICES. 
IFLQR = 3s FILE POSITION NUMBER OF IUTQR. 


t ORE CM E0E0M2 OYLOFLX 

(WAGNER .5 .5 1.0 2.0 


18. 2C 1 
18.3 ) 


■ OR WAGNER OPTION. 

AQR * .5COE+00* WAGENR FUNCTIONS Ao 

BQR * o500E+00t WAGENR FUNCTION* B. 

ALQR « . 100E+01 * WAGENR FUNCTION* ALPHA. 

BEQR a .200E+01 * WAGNER FUNCTION* BETA. 

( SQUIT 19.0 » 
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MATRIX EQUATIONS 


TITLE CHECK CASE 2, MODIFY DERIVATIVES FRDM CARDS » ADD SENSORS, OR TP 

INPUT MATRICES FREQUENCY l 


S F N S G R 


— 

— 

— matrix 

H2BAR 

DIMENSIONED I 


3 X 

A) 

ROW 

l 

0- 

0* 

0. 

0. 



ROW 

2 

0. 

0. 

0* 

0. 



ROW 

3 

.10006*01 

— .1 OOOE *02 . 3000E * 01 

0. 






M3BAR 

DIMENSIONED t 


3 X 

A) 





ROW 

l 

• 10006+01 

.12006+02 .8 0006+00 

0. 



ROW 

2 

• IQ00E+01 

-.60006*01 . 50006 * 00 

0. 



ROW 

3 

0. 

0. 

0. 

0. 



OUT 

P 

U T MATS 

ICES 

FREQUENCY i 







MATRIX 

MI 

DIMENSIONED I 

6X 

6) 



ROW 

1 

0. 

0. 

0. 

0. 


0. 

0. 

ROW 

2 

0* 

c. 

0. 

0. 


0. 

0. 

ROW 

3 

0. 

0. 

. 2000E* 03 

0. 


0. 

0. 

ROW 

A 

0. 

0. 

0. 

0. 


0. 

0. 

ROW 

5 

0* 

0. 

0. 

0. 


-.10006*01 

0. 

ROW 

6 

0. 

0. 

0. 

0. 


0. 

lOOOE+Oi 

— 

— 

MATRIX 

N 3 

DIMENSIONED ( 

6X 

61 



ROW 

1 

.IOOOE+OA 

0. 

0. 

0. 


0. 

0. 

ROW 

2 

0. 


AOOOE *05 0. 

0. 


0. 

0. 

ROW 

3 

0. 

0. 

.20006+02 

•5000E*QL 

0. 

0. 

ROW 

A 

0. 

0. 

0. 

.15006*02 

0. 

0. 

ROW 

5 

• 10006*01 


I200E+02 .80006*00 

Oo 


0. 

0. 

ROW 

6 

. lOCOE+Ci 

- • 

60C0E *01 .50006*00 

0. 


D. 

0. 

— 

— 

MATRIX 

FREQH DIMENSIONED l 

2 X 

n 




ROW I 0. 

ROW 2 .20006+01 
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INPUT 


MATRICES FREQUENCY 


1 


OUTPUT MATRICES 


FREQUENCY 


— — matrix 


M4 


DIMENSIONED I 


ex 


61 


ROM 

L 

Oo 

.34976+02 

' 60006+01 

.20016+02 

ROW 

2 

Oo 

-.49996*02 

.80006+01 

-.30006+02 

POM 

3 

Oo 

.40006+01 

.20006+01 

-.10006+01 

ROW 

4 

Oo 

o 1 000E+0 1 

10006+01 

.40006+01 

now 

5 

Oo 

Oo 

Do 

0. 

ROW 

6 

Oo 

0. 

0. 

Oo 

— 


- MATRIX 

M5 0IMENSI0N60 C 

6X 61 

ROW 

l 

.35016+01 

.10036+01 

.70006+01 

.30006+01 

ROW 

2 

-.50006*01 

-.30006+01 

.40006+01 

.20006+01 

ROW 

3 

—.40006*00 

.10006+01 

.20006+01 

.10006+01 

ROW 

4 

-.10006+00 

o 10006+01 

. 10006+ 01 

.20006+01 

ROW 

5 

0. 

0. 

0. 

0. 

ROW 

6 

0. 

0. 

0. 

0. 


0 . 
0 . 
Oo 
0. 
0 « 
Oo 
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OOOOOO oooooo 


input matrices FREQUENCY z 


— 

— 

- MATRIX 

M4 

DIMENSIONED f 

4X 4) 

ROW 

1 

0. 

— o3000E *-01 

« 7Q00E + 01 

.iOOOE+Ol 

ROW 

Z 

0* 

.6 0006+01 

• 90006+ 01 

• lOQOE+Ql 

ROW 

3 

0. 

. 5000E+01 

. 3000E+01 

-. 1000E+01 

ROW 

4 

0. 

•zocoe+oi 

LOOOE+ 01 

.3000E+0! 

— 

— 

- MATRIX 

M5 

DIMENSIONED C 

4X 4) 

ROW 

L 

.30CQE+Q0 

.60006+01 

. 600QE+ 01 

.4000E+01 

ROW 

2 

— .60CQE+00 

.5000E+01 

. 3000E+01 

.3000E+01 

ROW 

3 

-.SOOOE+OO 

-2O0OE+OI 

. 1000E+01 

.2000E+01 

ROW 

4 

-.2000E+00 

• 20006 + 01 

0. 

•2000E+01 

— 

— 

- MATRIX 

C3 I 

DIMENSIONED I 

8 X 1) 

ROW 

i 

.30006+00 




ROW 

2 

.iogoe+oo 




ROW 

3 

- » bO OCE +00 




ROW 

4 

.2000E+00 




ROW 

5 

-,5000E+00 




ROW 

6 

. LOQOE+OO 




ROW 

7 

-.2000E+00 




ROW 

8 

. 3000E + 00 




OUT 

P U 

T M A T R 

ICES FREQUENCY 2 



ROW 

I 

MATRIX 

0. 

H4 1 

.34576+02 

DIMENSIONED { 
. 7000E+ 01 

6 X 6) 

.200 1E + 02 

0. 

0. 

ROW 

2 

0. 

-.4999E+02 

= 9000E + 01 

— . 3000 E *02 

0. 

0. 

ROW 

3 

0* 

• 5Q00E + 0L 

. 3000E + OL 

-« 1000E+01 

Oo 

0. 

ROW 

4 

0. 

.20006+01 

loooe + oi 

.30006+01 

0. 

0. 

ROW 

5 

0. 

0. 

Oo 

0. 

0. 

0. 

ROW 

6 

0. 

0. 

0. 

0. 

0. 

0. 

ROW 

I 

MATRIX 

.35016+01 

M 5 DIMENSIONED I 

• 1 0C3E +0 1 e&OOQE+Ol 

6X 6) 

.4000E+01 

0. 

0. 

ROW 

2 

-.50006*01 

-.30006+01 

o 3000E+01 

. 3000E+01 

0. 

0. 

ROW 

3 

- .50006+00 

.20006+01 

. LOOOE+OI 

.2000E+01 

0. 

0. 

ROW 

4 

-.20006*00 

. 2000E + 01 

Oo 

.20006+01 

0. 

0. 

ROW 

5 

0. 

0. 

0. 

0. 

0. 

0. 

ROW 

b 

0. 

0. 

0. 

0. 

Oo 

0. 
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LOADS MATRIX EQUATIONS 

TITLE CHECK CASE 2* MODIFY DERIVATIVES FROM CARDS* ADD SENSORS. QR TP 

INPUT MATRICES FREQUENCY l 

OUTPUT MATRICES FREQUENCY l 


81 



I N P U 
0 U T P 

ROW 1 
ROW 2 


T MATRICES FREQUENCY l 
UT MATRICES FREQUENCY l 

— MATRIX M4BAR DIMENSIONED ( 

• 6000E+01 •4000E+01 .3000E+01 

.60006+0 l .3000E+0 1 .5000E+01 


2 X 61 

0 . 0 . 
•2000E+01 0. 
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INPUT MATRICES FREQUENCY 2 


— — MATRIX. M4BAR DIMENSIONED I 2X 41 

ROW I •7000E+0 1 •3000E+01 .2000E«-0l -.IOOOE+Ol 

ROW 2 •50C0E+01 .2000E+01 . 4000E+01 . IOOOE+Ol 

— — — t— p- MATRIX M5BAR DIMENSIONED ( 2X 4) 

ROW l ^.30C0E*Cl 0. -« 2000E+ 01 .4000E+01 

ROW 2 .2COOE+01 0. -.3000E+01 -.IOOOE+Ol 

MATRIX C3BAR DIMENSIONED ( 4X II 

ROW I -.30C0E+0I 

ROW 2 • lOOOE+CO 

ROW 3 . 20C0E+C I 

ROW 4 « 20C0E+00 

OUTPUT MATRICES FREQUENCY 2 


— - MATRIX M4BAR DIMENSIONED I 2X 6) 

ROW 1 • 1400E+02 .60Q0E+01 .4000E+01 -.2000E+01 0. 

ROW 2 .1000E+02 .4000E+01 . 8000E + 01 .2000E+01 0. 
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QR MATRICES GENERATED 
ROOT OPTION 


ROW 

ROW 

ROW 

ROW 

ROW 

ROW 


ROW 

ROW 

ROW 

ROW 

ROW 

ROW 


ROW 

ROW 

ROW 

ROW 

ROW 

ROW 


TITLE CHECK CASE 2, MODIFY DERIVATIVES FROM CAROS, ADD SENSORS* QR TP 


MATRIX 


S**2 


DIMENSIONED < 


6 X 


6) 


• 10 OOE+G4 

0. 

0. 

0. 

. 1000E+01 


0. 

•4000E+05 

0. 

0. 

•1200E+02 


• lOOOE+Ol — .6000E+0 1 


0. 

0. 

.20006+02 

0. 

. 8000E+ 00 
• 50006+ 00 


0. 

0. 

•5000E+0! 
• 1500E+02 
0. 

0. 


0. 

0. 

0. 

0. 

0. 

0. 


MATRIX S**l 


DIMENSIONED ( 6X 6) 


1 

2 

3 

4 

5 

6 


.35016+01 
— .5000E+01 
-.4000E+00 
-.10006+00 
0. 

0. 


.10036+01 

-.30006+01 

.10006+01 

.10006+01 

0. 

0. 


. 70006+01 
• 4000E + 01 
. 2000E + 01 
.10006+01 
0. 

0. 


.30006+01 0. 
•2000E+01 0. 
• 1000E+01 0. 
•2000E+01 0. 

0 . 0 . 

0 . 0 . 


0. 

0. 

0. 

0. 

0. 

0. 


MATRIX S*40 


DIMENSIONED ( 6X 6) 


L 0. 

2 0 . 

3 0. 

4 0. 

5 0. 

6 0 . 


.34976+02 

-.4999E+02 

•4000E+01 

•1000E+01 

0. 

0. 


.60006+01 

• 8000E + 01 

• 2020E + 03 
1000E+01 

0. 

0. 


•2001 E+02 
3000E+02 
-.10006+01 
.40006+01 

0. 

0. 


0. 

0. 

0. 

0. 

-. 1000E+01 

0. 


0. 

0. 

0. 

0. 

0. 

-• lOOOE+Ol 
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OR MATRICES GENERATED 
TIME SOLUTION OPTION 


TITLE CHECK CASE 2, MODIFY DERIVATIVES FROM CARDS, ADD SENSORS, QR TP 
MATRIX S**2 DIMENSIONED ( BX 81 


ROW 

1 

• 10G0E+C4 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

ROW 

2 

0. 

.40ooe*o5 

0. 

0. 

0. 

0. 

0. 

0. 

ROW 

3 

0. 

0. 

. 2000E* 02 

. 5000E+01 

0. 

0. 

0. 

0. 

ROW 

4 

0. 

0. 

0. 

• 1500E+02 

0. 

0. 

0. 

0. 

ROW 

5 

• 1000E*,01 

.12008*02 

• 8000E* 00 

0. 

0. 

0. 

0. 

0. 

ROW 

6 

.10006*01 

- .60006*01 

. 5000E+00 

0. 

0. 

0. 

0. 

0. 

ROW 

7 

.7000E+01 

. 8C006+01 

• 4000E* 01 

• 1000E+01 

0. 

0. 

0. 

0. 

ROW 

8 

•60C06+01 

.30006+01 

lOOOE+Ol 

.2000E+01 

0. 

0. 

0. 

0. 

ROW 

I 

MATRIX 

.35C1E+01 

S**l DIMENSIONED ( 

. 1003E *01 . 7000E+ 01 

8 X 81 

.30006+01 

0. 

0. 

0. 

0. 

ROW 

2 

-.5000E+01 

-•3000E + 01 

. 4Q0QE+ 01 

.2000E+01 

0. 

0. 

0. 

0. 

ROM 

3 

-.4000E+00 

. lOOOE + Ol 

. 2000E+01 

. 1000E+01 

0. 

0. 

0. 

0. 

ROW 

4 

-.lOOOE+OO 

. lOOOE + Ol 

. IOOOE* 01 

.20006+01 

0. 

0. 

0. 

0. 

ROW 

5 

0. 

0. 

0. 

-0. 

0. 

0. 

0.. 

0. 

ROW 

6 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

ROW 

7 

-.2000E+01 

.10006*01 

IOOOE* 01 

• 5000E+01 

0. 

0. 

0. 

0. 

ROW 

8 

.3000E+01 

- .lOOOE+Ol 

— • 2000E* 01 

-.20006*01 

0. 

0. 

0. 

0. 

ROW 

l 

MATRIX 

0. 

S*40 DIMENSIONED ( 

•3497E+02 • 6000E+ 01 

8X 81 

.2001E+02 

0. 

0. 

0. 

0. 

RON 

2 

0. 

— .4999E *02 

. 8000E * 01 

-.3000E+02 

0. 

0. 

0. 

0. 

ROW 

3 

0. 

.4000E *01 

. 20206+ 03 

-. 1000E+01 

0. 

0. 

0. 

0. 

ROW 

4 

0. 

. 1000E*0 1 

IOOOE* 01 

•4000E+01 

0. 

0. 

0. 

0. 

ROW 

5 

0. 

0. 

0. 

0. 

1000E+0 1 

0. 

0. 

0. 

ROW 

6 

0. 

0. 

0. 

0. 

0. 

-. 1000E+01 

0. 

0. 

ROW 

7 

• 80 00 E <-01 

.400QE *0 1 

. 3000E+01 

0. 

0. 

0. 

— • 10006+01 

0. 

ROM 

8 

• 60 COE + 01 

«3000E*01 

.50006+01 

.20006*01 

0. 

0. 

0. 

-. L000E+01 


MATRIX VECTOR DIMENSIONED I 8X 1) 


RON 1 — • 350 1E+01 

ROW 2 .50006*01 

ROW 3 -.4000E+00 

ROW 4 -.10006*00 

ROW 5 0. 

ROW 6 0. 

ROW 7 - .20QQE+QI 

ROW 8 . 3000E + 01 
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OR MATRICES G E N £ R A T E D < 
WAGNER FUNCTION OPTION 


^ _ 


T! 

MATRIX 

: TLE CHECK CASE 2* MODIFY DERIVATIVES 

S**4 DIMENSIONED I 6X 6) 

FROM CARDS* ADO SENSORS* QR TP 

ROW 

1 

.1000E+04 

0. 

0. 

0. 

0. 

0* 

ROW 

2 

0. 

•4000E+05 

0. 

0. 

0. 

0. 

ROW 

3 

0. 

0. 

• 2000E+-02 

• 50Q06+01 

0. 

0. 

ROW 

4 

0. 

0. 

0. 

• 1 500E*02 

0. 

0. 

ROW 

5 

0. 

0. 

0. 

0. 

0. 

0. 

ROW 

6 

0. 

0. 

0. 

0. 

0. 

0. 

ROW 

1 

MATRIX 

•3000E+C4 

S**3 DIMENSIONED C 

0. 0. 

6 X 6) 

0. 

0. 

0. 

ROW 

2 

0. 

• I200E+06 

0. 

0. 

0. 

0. 

ROW 

3 

0. 

0. 

• 6000E*- 02 

•1500E+02 

0. 

0. 

ROW 

4 

0. 

0. 

0. 

.4500E*02 

0. 

0. 

ROW 

5 

0. 

0. 

0. 

0. 

0. 

0. 

ROW 

6 

0. 

0. 

C. 

0, 

0. 

0. 

ROW 

1 

MATRIX 

„2005E*04 

S**2 DIMENSIONED I 

» L 504E «-01 • I050E+ 02 

6 X 6) 

•4500E+01 

0- 

0. 

ROW 

2 

-.7500E+0I 

-8000E+05 

• 6000E+ 01 

• 3000E+01 

0. 

0. 

ROW 

3 

-,60COE*00 

• 1 500E *01 

• 2430E + 03 

• 1 1 50E + 02 

0. 

0. 

ROW 

4 

— • 1500E+00 

•1500E+01 

. I500E+ 01 

.3300E+02 

0. 

0. 

ROW 

5 

• 10C0E + G1 

.1200E+02 

. 8000E+00 

0. 

0. 

0. 

ROW 

6 

-1000E+01 

-•6000E+01 

. 5000E+00 

0. 

0. 

0. 

ROW 

I 

MATRIX 

.7002E+01 

S**I D 

• 544 7E *02 

I ME NS I ONE 0 C 
.2300E+02 

6 X 61 

.3601E+02 

0. 

0. 

ROW 

2 

-.1000E*02 

-.8099E«-02 

•2000E+02 

-.4100E+02 

0, 

0. 

ROW 

3 

-.80C0E+00 

. 8000E «-0 I 

.60 70E + 03 

.50006*00 

0. 

0. 

ROW 

4 

-•2000E+00 

.3500E + 01 

. 5000E«-00 

. 10006*02 

0. 

0. 

ROW 

5 

0. 

0. 

0. 

0. 

0. 

0. 

ROW 

6 

0. 

0. 

0. 

0. 

0. 

0. 

ROW 

I 

MATRIX 

0. 

S440 DIMENSIONED ( 

.6995E*02 . 1200E + 02 

6 X 6) 

.40026*02 

0. 

0. 

ROW 

2 

0. 

— ,99986 *02 

. 1600E* 02 

-.60006*02 

0. 

0. 

ROW 

3 

0. 

.SOOOE^OI 

• 4040E*-03 

—.20006*01 

0. 

0. 

ROW 

4 

0. 

. 2000E *-01 

-•2000E4-01 

. 80006*01 

0. 

0. 

ROW 

5 

0. 

0. 

0. 

0. 

10006*0 l 0. 

ROW 

6 

0. 

C. 

0. 

0. 

0. 

-. 1000E*01 
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•*< 

• 

• 

PROGRAM 

L219A2 

VERSION JUNE 29,77 IS FINISHED- 

** 

• 

DATE OF 

RUN IS 

77/11/09. 


• 

TIME OF 

RUN IS 

09-46.22- 

** 


I *** * ****** * ********** ******************************* 


PROGRAM L219A2 VERSION JUNE 29,77 NOW RUNNING. 
THE PROGRAM IS PART OF THE DYLOFLX SYSTEM 
DEVELOPED FOR NASA UNDER CONTRACT NAS1-13918. 
DATE OF RUN IS 77/U/09. 

TIME OF RUN IS 09.46.23. 


******* ******** ****** ******** ***** ************ ******* 

( *EOMOD TEST CASE 3 

ITITLE CHECK CASE 3, MODIFY FOM AND LOADS 
(SIZE 4 2 


1.0 


3.0 


> 

I 

I 


P R 0 B L 
NDOF 
NP AN 

NFREOM * 
l OUTPUT TEST4 


EM SIZE 

4. TOTAL NUMBER OF DEGREES OF FREEDOM. 
0, TOTAL NUMBER OF PANELS. 

2, NUMBER OF FREQUENCIES. 

I DUMMY 1 


4.1 » 


( PRINT 


OUTPUT 
IUTEOM = TEST4 
IFLEOM « 

IUTLOD » DUMMY 
IFLLOO * 

OUTPUT 


TAPES 

• TAPE NAME FOR MODIFIED EOM MATRICES. 

1, FILE POSITION NUMBER OF IUTEOM. 

, TAPE NAME FOR MODIFIED LOADS MATRICES. 
1, FILE POSITION NUMBER OF IUTLOD. 
MATRICES ALL 


4.3 I 


OUTPUT MATRICES PRINT OPTION. 

IUTPR = —999, IF IUTPR = -999, PRINT ALL OUTPUT MATRICES. 

= 0, NO OUTPUT MATRICES PRINTED. 

= N, MATRICES OF NTH FREQUENCY ONLY PRINTEO. 

= 999, PRINT ONLY MODIFIED MATRICES. 

I SYMMETRIC 5.0 ) 

ISEOM EOM TAP 2DYL0FLX 6.0 ) 

EQUATIONS OF MOTION 


INEOM = EOMTAP , EOM INPUT TAPE NAME. 

INEOMF = 2. FILE POSITION NUMBER OF INEQM. 


(SCALE M 4 


2.0 



8. 0 

) 

(BODY AXIS 1 

2 3 

1.0 

10.0 


13.0 

) 

BODY AXIS 

SYMMETRIC 






I COL X « 

W 

X- COLUMN 

OF ORIGINAL 

AXIS. 



ICOLZ * 

2, 

Z- COLUMN 

OF ORIGINAL 

AXIS. 



ICOLT = 

3, 

THE TA-COLUMN 

OF ORIGINAL 

AXIS. 



ALPHA! * 

.L75E-0I, 

ALPHAl, ANGLE 

OF ATTACK. 




BODYVT » 

. 100E*02, 

VT, VELOCITY 

- TRUE AIR 

SPEED. 



($LOADS LOOTAP 

2 

20YLOFLX 



14.0 

1 


LOADS EQUATION 


INLOO * LODTAP , LOAOS INPUT TAPE NAME. 

INLODF = 2. FILE POSITION NUMBER OF INLOD. 

NLDOU = 2, NUMBER OF OUTPUT LOADS. 


(REPLACE M36AR 
(REPLACE M4BAR 
( INCREMENT M JBAR 
( iQUIT 


2 


I 4 10. 0 

1 4 5. 0 

2 4 3.0 


16.0- 1 » 
16.0-2 I 

17.0 I 

19.0 I 
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EOM MATRIX 


EQUAT IONS 





TITLE 

CHECK CASE 3. MODIFY 

EOM 

AND LOADS 

OUT 

R 

U T HAT 

R 

I C E 

S FREQUENCY 1 





— 

MATRIX 


Ml 

DIMENSIONED ( 


4X 

41 

ROW 

1 

O. 


0. 

0. 

0. 



ROW 

2 

0. 


0. 

0. 

0. 



ROW 

3 

0. 


0. 

. 1000E+03 

0. 



ROW 

4 

0. 


0. 

0. 

0. 



— 

— 

— MATRIX 


M2 

DIMENSIONED ( 


4X 

4 1 

ROW 

1 

0. 


0. 

. I745E»03 

0. 



ROW 

2 

0. 


0. 

— • 4000E+ 06 

0. 



ROW 

3 

0. 


0. 

0. 

0. 



ROW 

4 

0. 


0. 

0. 

0. 



— 

— 

MATRIX 


M3 

DIMENSIONED I 


4 X 

4) 

ROW 

1 

• 1000E+C4 

0. 

0. 

0. 



ROW 

2 

0. 


.4000E+05 0. 

0. 



ROW 

3 

0. 


0. 

• I000E+02 

0. 



ROW 

A 

0. 


0. 

0. 

• 

1 500E+02 


— 

— MATRIX 


FREQM 

DIMENSIONED I 


2 X 

l) 

ROW 

1 

0. 







ROW 

2 

•2000E+01 
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OUTPUT MATRICES FREQUENCY i 


— 

— 

— MATRIX 

M4 1 

DIMENSIONED 

ROW 

1 

O. 

-.40006*01 

-.57976402 

ROW 

2 

0 - 

.1 0006*02 

-.4409E402 

ROW 

3 

0. 

•80C0E401 

-.607064 01 

ROW 

4 

0. 

•2000E 401 

-.12026402 

— 

— 

-- MATRIX 

M5 DIMENSIONED 

ROW 

1 

.20006400 

.70006401 

.70006401 

ROW 

2 

-.5000E400 

•60006401 

.40006401 

ROW 

3 

-.40006400 

•1000E401 

• 20006 + 01 

ROW 

4 

-.lOOOE+OO 

•10006401 

.10006+01 

— 

— 

— MATRIX 

C3 DIMENSIONED 

ROW 

1 

•200CE+00 



ROW 

2 

0 . 



ROW 

3 

-.50COE400 



ROW 

4 

0 . 



ROW 

5 

-.40C0E4C0 



ROW 

6 

0 . 



ROW 

7 

-.1000E400 



ROW 

8 

0. 




4 X 4) 

•40006+01 
.2000E+0 l 
-• 2000E+01 
• 8000E+01 

4 X 4) 

•30006*01 
.2000E+01 
.lOOOE+OI 
• 2000E+01 

8 X 1) 



OUTPUT MATRICES 


FREQUENCY 


2 


— 

— 

MATRIX 

M4 

DIMENSIONED l 

4 X 41 

ROW 

1 

O. 

-.30006+01 

-.5295E + 02 

• 1000E+01 

ROW 

2 

0. 

.60006+01 

-.41106+02 

.10006+01 

ROW 

3 

0. 

*50006+01 

-• 1709E+02 

— . 1000E+01 

ROW 

4 

0. 

•2000E+01 

— • 21 03E+02 

.3000E+01 

— 

— 

MATRIX 

M5 

DIMENSIONED 1 

4 X 4) 

ROW 

1 

• 3000E+00 

•60006+01 

■ 60006 + 01 

•4000E+01 

ROW 

2 

-•60006+00 

• 5000E + 01 

• 3000E+ 01 

• 3000E+01 

ROW 

3 

-•50006+00 

.20006+01 

• 1000E+01 

.20006+01 

ROW 

4 

-•2000E+00 

•20006+01 

0. 

.20006+01 

ROW 

ROW 

ROW 

ROW 

ROW 

ROW 

ROW 

ROW 

1 

Z 

3 

4 

5 

6 
1 
e 

MATRIX 

•30006+00 
• 1000E+00 
-.6QCOE+CO 
•2000E+CO 
-.5000E+CO 
• lOOOE+CO 
-•20006+00 
•30006+00 

C3 

DIMENSIONED I 

ax i) 



LOADS MATRIX 


EQUATI ONS 


TITLE CHECK CASE 3 t MOOIFY EON AND LOADS 


OUT 

P U 

T M A T R 

ICES 

FREQUENCY I 





— 

MATRIX 

MIBAR 

DIMENSIONED ( 


2 X 

A) 

ROM 

1 

O. 

0 . 

0 , 

0 . 



ROM 

2 

0 * 

0 . 

0 . 

0 . 



— 

MATRIX 

M2BAR 

DIMENSIONED 1 


2 X 

41 

ROM 

I 

0 . 

0 . 

— • 7878E* 02 

0 . 



ROM 

2 

0 . 

0 . 

“•2895E*'02 

0 . 



— 

— 

MATRIX 

M3BAR 

DIMENSIONED ( 


2 X 

41 

ROM 

1 

•7000E+0 l 

.8000E+01 . 4000E+01 

• 

1000E+02 

ROM 

2 

.6C0OE+01 

.3000E+01 -• 1000E+0I 

. 

5000E+01 
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OUTPUT MATRICES FREQUENCY 1 


MATRIX MABAR DIMENSIONED t 2X A1 

ROW l . 8000E+0I .AOOOE+Ol 73A9E*0l 0. 

ROW 2 •6000E+0 l •3000E+0! . 1552E+02 .2000E+01 

MATRIX M5BAR DIMENSIONED I 2X AJ 

ROW l -’•20GQE+OI .1000E + 01 -.lOQOE+Oi -5000E+01 

ROW 2 .3000E+01 -.1000E*0l -.2000E+0! -.2000E+01 

MATRIX C3BAR DIMENSIONED ( 4X I) 

ROW 1 -•2000E+C1 

ROW 2 0. 

ROW 3 • 3000E+0 1 

ROW A 0. 
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OUTPUT MATRICES 


FREQUENCY 


2 


MATRIX 


MABAR 


DIMENSIONED I 


2 X 


A ) 


ROW 

ROW 


ROW 

ROW 


• 7000E+01 
• 5000E+0I 

MATRIX 


.3000E + OI 
.2000E + 01 


M5BAR 


-•3000E+CI 0. 

. 2000E+C l 0. 


• 1AT6E*-0I 
. A349E* 01 

DIMENSIONED C 


■ 5000E+01 
« 1000E+01 


2 X 


Al 


-.2000E*01 .A000E+01 

-.3000E4-0L -• LOOOE+OI 


MATRIX 


C3BAR 


DIMENSIONED I 


AX 


l) 


ROW 1 
ROW 2 
ROW 3 
ROW A 


-.3000E+01 
• I0CCE+00 
•2000E+0! 
.2000E+00 
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PROGRAM 121**2 VERSION JUNE 29*77 IS FINISHED. 
DATE OF RUN IS 77/11/09. 

TINE OF RUN IS 09.46.25. 


I*************************************************** 


************ ****************** ******** ************** 


PROGRAM L219A2 VERSION JUNE 29,77 NOW RUNNING. • 

THE PROGRAM IS PART OF THE DYLOFLX SYSTEM * 

DEVELOPED FOR NASA UNDER CONTRACT NASi-13918. * 

DATE OF RUN IS 77/11/09. * 

TIME OF RUN IS 09.46.26. * 

***************************************************** 

(SEQMOD TEST CASE 5 1.0 > 

(TITLE CHECK CASE 5, MODIFY DERIVATIVES FROM CARDS AND ADD SAS AND SENSORS » 

(SIZE 7 2 3.0 I 


P R 0 E L 
NDOF 

NP AN = 
NFREQM = 

( OUT FUT EOEOM 


EM SIZE 

7* TOTAL NUMBER OF DEGREES OF FREEDOM. 
0, TOTAL NUMBER OF PANELS. 

2, NUMBER OF FREQUENCIES. 

1 EQLOD 1 


4-1 | 


(PRINT 


OUTPUT TAPE 
IUTECM « EQEOM , 

IFLEOM * 1* 

IUTLOD * EQLOD * 

IFLLOD * 10* 

INPUT MATRIX 


S 

TAPE NAME FOR MOOIFIED EOM MATRICES. 
FILE POSITION NUMBER OF IUTEOM. 

TAPE NAME FOR MODIFIED LOADS MATRICES. 
FILE POSITION NUMBER OF IUTLOD. 

ALL 


4.2 ) 


INPUT MATRICES PRINT OPTION. 

INPR * -999, IF INPR = -999, PRINT ALL INPUT MATPICES. 

= 0, NO INPUT MATRICES PRINTED. 

= N, MATRICES OF NTH FREQUENCY ONLY PRINTED. 

(PRINT OUTPUT MATRIX ALL 4.3 ) 


OUTPUT MATRICES 
IUTPR * -999, IF 


( SYMMETRIC 
( SEOM EOMTAP 


PRINT OPTION. 

IUTPR » -999, PRINT ALL OUTPUT MATRICES. 

= 0, NO OUTPUT MATRICES PRINTED. 

= N, MATRICES OF NTH FREQUENCY ONLY PRINTED. 

= 999, PRINT ONLY MODIFIED MATRICES. 

5.0 

1 1 0 1 l 1 1 0 0 6.0 


I 

> 


EQUATIONS OF 


M 0 T I ON 


IN EOM * EOMTAP , EOM INPUT TAPE NAME. 

INEOMF * 1, FILE POSITION NUMBER OF INEOM. 4 . 

(DERIVATIVE FROM CARDS 1 0 2.0 10.0 7.1 I 


DERIVATIVES 

I VOL 

NCS 

I NOUN « 
QUEBAR * . 

VT «= . 


FOR SYMMETRIC ANALYSIS. 

CARO, INPUT VOLUME 

1, NUMBER OF CONTROL SURFACE 
0, UNSTEADY DERIVATIVE INDICATOR 
200E*01, DYNAMIC PRESSURE 
100E+02, VELOCITY (TRUE AIR SPEED) 


COLUMN NUMBERS OF RIGID BODY FREEDOMS 
IXCOL * 0, COLUMN OF X FREEDOMS 

IZCOL - 1, COLUMN OF Z FREEDOMS 

ITCOL * 2, COLUMN OF THETA FREEDOMS 

COLUMN OF DELTA CONTROL SURFACE FREEDOMS FOLLOW 
4 


CONSTANTS ASSOCIATED WITH DERIVATIVES 

XREF = 0. t X-COORDINATE OF MOMENT REFERENCE 
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ZREF 

0. t 

Z-COORDINATE OF MOMENT REFERENCE 

ALPHA 1 * 

• I 75E-02t 

IG ANGLE OF ATTACK 

SW - 

• 2 50E+C L v 

WING REFERENCE AREA 

C BAR * 

• 1 OOE + O 1» 

REFERENCE CHORD 

CLIP 

•3CCE+01, 

RIGID STEADY STATE DERIVATIVE 

CL IE 

•20CE+00t 

ELASTIC STEADY STATE DERIVATIVE 

STEADY STATE DERIVATIVES FROM CARD 

CLU 

.4C0E+0U 

C-L-UHA T-PIGID 

CDU 

. 1 COE+ 0 1 y 

C-O-UHAT-RIGID 

CMUREF 

•lOOE+Olt 

C-M-UHAT-RI GI D-REF 

CLA 

• 700E+0 l t 

C—L— ALPHA— RIGID 

CDA 

. 1006+01* 

C-D- ALPHA-RIGID 

CM ARE F 

• 100E+02* 

C-M-ALPHA-R1G1D-PEF 

CLUE 

• 1 006+0 1 • 

C-L-UHA T-ELASTIC 

CDUE 

• 1 OOE + OOt 

C-O-UHAT-ELASTIC 

CLAE 

-• lOOE+Olt 

C-L-ALPHA-ELASTIC 

COAE 

• 100E+0 It 

C-D- ALPHA- ELASTIC 

CLO 

.200E+01» 

C-L-QUEHAT-rRIGID 

CDQ 

. 300E+0 l» 

C— D— QUE HA T— R I G I 0 

CMOREF 

• 7 OOE + O It 

c-n-quehat-rigio-ref 

CLQE 

• lOOE+OOt 

C-L-OHAT-ELASTIC 

CDQE 

-. 100E«-01, 

C-D-OHA T-ELASTIC 


CONTROL SURFACE DERIVATIVES 
CLO (C-L-OELTA-R IGIDI FOLLOWS 
•400E+01 

COO (C-O-OELTA-R IGIDI FOLLOWS 
• 100E+01 

CMOREF ( Q-H-OELTA-R BF I FOLLOWS 
•60CE+0 I 

CLDE (C-L-DELTA-ELASTIC) FOLLOWS 
-.200E+01 

CDOE ( C— O-DEL T A- EL A STIC I FOLLOWS 
-.100E+01 

(SENSOR 3 0 11 9.1 ) 


INSEN - L00TP2 » SENSORS INPUT TAPE NAME. 

INSENF « It FILE POSITION NUMBER OF INSEN. 

NLOSEN « 3» NUMBER OF LOADS ON SENSOR TAPE. 

MATRIX IN OUT 

( M2BAR 3 5 

(M3BAR 1 6 

( ENO SENSOR 
(SAS 


9.2- 1 t 

9.2- 2 ) 

9.3 I 

10.1 ) 


ROW 

COL Mi 

M2 

M3 



( 7 

7 1.0 



10.2-1 


( 7 

5 

-1.0 


10.2-2 


( 7 

6 


2.0 

10.2-3 


( 4 

7 1.0 

2.0 

3.0 

10.2-4 


(END SAS 




10.3 



( SOU IT ) 


97 











E 

0 M 

MATRIX E Q 

U A T I 

0 N S 



TITLE 

CHECH CASE 5t MODIFY DERIVATIVES FROM CARDS ANO ADD SAS AMO SENSORS 

I N P 

U 

T M A T R I 

</> 

Ul 

o 

FREQUENCY 1 



— 

MATRIX 

Ml 

DIMENSIONED ( 

AX 

At 

SOW 

1 

0. 

0* 

0. 

0. 


sow 

2 

0. 

0. 

0. 

0. 


ROW 

a 

0. 

0. 

. 1000E+03 

0. 


ROW 

A 

0. 

0. 

0. 

0. 


— 

— 

MATRIX 

M3 

DIMENSIONED ( 

AX 

Al 

ROW 

1 

• 10Q0E+CA 

0. 

0. 

0. 


ROW 

2 

0. 

.A000E+05 0. 

0. 


ROW 

3 

0. 

0. 

« 1000E* 02 

0. 


ROW 

A 

0. 

0. 

0. 

•1500E+02 



S E N S 0 

R 




— 

MATRIX 

M2BAR 

DIMENSIONED A 

3X 

Al 

ROW 

1 

0. 

0. 

0. 

0. 


ROW 

2 

0. 

0. 

0. 

0. 


ROW 

3 

• ICC0E+C1 

- . 1000E+02 ■ 300QE+ 01 

0. 


— 

— 

■ — MATRIX 

M3BAR 

DIMENSIONED I 

3X 

Al 

ROW 

I 

.lOOOE+Oi 

• 1200E+Q2 . 80 00E* 00 

0. 


ROW 

2 

•lOOOE+Oi 

-.600QE + 01 * 5000E + 00 

0 . 


ROW 

3 

0. 

0. 

0 . 

0 . 


OUT 

P 

U T M A T R 

I C E 

S FREQUENCY l 




— 


MATRIX 

MI 

DIMENSIONED ( 


7X 

71 




ROW 

1 

0. 

0. 

0. 

0. 



0. 

0. 

0. 

ROW 

2 

0. 

0. 

0. 

0. 



0. 

0. 

0. 

ROW 

3 

0. 

0. 

. 1000E* 03 

0. 



0. 

0. 

0. 

ROW 

A 

0. 

0. 

0. 

0. 



0. 

0. 

.10006*01 

ROW 

5 

0. 

0. 

0. 

0. 



1000 E*0 1 

0. 

0. 

ROW 

6 

0. 

0. 

Q. 

0. 



0. 

-.10006*01 

0. 

ROW 

7 

0. 

0. 

0. 

0. 



0. 

0. 

a lOOOEwOl 

— 

— 

MATRIX 

M2 

DIMENSIONED ( 


TX 

71 




ROW 

1 

0. 

0. 

0. 

0. 



0. 

0. 

0. 

ROW 

2 

0. 

0. 

0. 

0. 



0. 

0. 

Oa 

ROW 

3 

0. 

0. 

0. 

0. 



0. 

0. 

0. 

ROW 

A 

0. 

0. 

0. 

0. 



0. 

0. 

. 2000E*01 

ROW 

5 

• 1000E + 01 

■“ « 

10006*02 . 3000E+ 01 

0. 



0. 

0. 

0. 

ROW 

6 

0. 

0. 

0. 

0. 



0. 

0. 

0. 

ROW 

7 

0. 

0. 

0. 

0. 



1Q00E+0 1 

0. 

0. 

— 

— 

MATRIX 

M3 

DIMENSIONED ( 


IK 

71 




ROW 

1 

. lOOOt+CA 

C. 

0. 

0. 



0. 

0. 

0. 

ROW 

2 

0. 

. 

AC00E+05 0. 

0. 



0. 

0. 

0. 
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ROM 

3 

0 . 

0 . 

• IOOOE+ 02 

0 . 


0 . 

0 . 

0 . 

ROM 

4 

0 . 

0 . 

0 . 

.15 00E+02 

0 . 

0 . 

.30006*01 

ROM 

5 

0 . 

0 . 

0 . 

0 . 


0 . 

0 . 

0 . 

ROM 

6 

» 1000E+01 

• 1200E*02 

. 8000E*00 

0 . 


0 . 

0 . 

0 . 

ROM 

7 

0 . 

- MATRIX 

0 . 

FREQM 

0 . 

dimensioned I 

0 . 

2 X 

l) 

0 . 

• 2000 E*0 1 

0 . 


ROM 1 O. 

Aon 2 • zoooe + o i 
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I N P 

U 

T 

M A T R I 

C E S FREQUENCY 1 


— 

— 

- MATRIX 

N4 

DIMENSIONED C 

4 X 4) 

ROW 

1 


0. 

- .2000E *01 

.6000E+01 

•2000E+01 

ROW 

2 


0. 

.5000E«-0l 

. 8000E ♦ 01 

* lOOOE+Ol 

P0W 

3 


0. 

•4000E+01 

• 2000E + 01 

-• 1000E+0I 

ROW 

4 


0. 

•l OOOE+OI 
M5 1 

1000E«-01 
DIMENSIONED I 

• 4000E* , 01 
4 X 41 




i n i a 

ROW 

L 


• 2000E+00 

» 7000E + 0 1 

. 7000E+01 

.3000E+01 

ROW 

2 


-• 5000E+00 

•6000E+0 1 

.4000E + 0! 

•2000E*0l 

ROW 

3 


-.4000E+00 

-1000E+01 

. 2000E*01 

• lOOOE+Ol 

ROW 

4 


IO0OE*OO 

• 1000E«-0l 

. 1000E + 01 

. 2000E+0 l 





C3 DIMENSIONED I 

8 X n 




” n A 1 K l A 

ROW 

1 


• 20CCE+00 




ROW 

2 


0. 




ROW 

3 


-,5000E*00 




ROW 

4 


0. 




ROW 

5 


-•4000E+00 




ROW 

6 


0* 




ROW 

7 


— • lOCOE+OO 




ROW 

8 


0. 




OUT 

P 

U 

T M A T R 

ICES FREQUENCY 1 



ROW 

1 

MATRIX 

0. 

M4 DIMENSIONED 1 

.3497E«-02 . 6000E ♦ 01 

7 X 71 

.2001E+02 

ROW 

2 

0. 

-.4999E*02 * 8000E* 01 

~.3000E*02 

ROW 

3 

0. 

•4000E*’OL • 2000E + 01 

-.lOOOE*Ol 

ROW 

4 

0. 

• lOCOE+Ol lOOOE+OI 

•4000E+01 

ROW 

5 

0. 

0. 0. 

0. 

ROW 

6 

0. 

0. 0. 

0. 

ROW 

7 

0. 

0. 0. 

0. 

— 

— 

MATRIX 

MS DIMENSIONED I 

7 X 71 

ROW 

1 

■ 350IE+0 1 

-5013E+00 .7000E+01 

• 3000E+01 

ROW 

2 

— • 5000 E*0 1 

— ■ l 750E* 01 « 4000E+ 01 

.2000E+01 

ROW 

3 

-.4000E*00 

. 1000E«-0 1 . 2 OOOE* Oi 

a IOOOE*OL 

ROW 

4 

— •1000E+00 

alOOOE+Ol a 1000E+01 

# 2000E*0l 

ROW 

5 

0. 

0. 0. 

0. 

ROW 

6 

0. 

0. 0. 

0. 

ROW 

7 

0. 

MATRIX 

0. 0. 

C 3 DIMENSIONED f 

0. 

1 4 X 1) 



ROW 

1 

- • 3501E+01 



ROW 

2 

0. 



ROW 

3 

•50O0E+0I 



ROW 

4 

0. 



ROW 

5 

— .4000E+00 



ROW 

6 

0. 



ROW 

7 

- • 1000E ♦OO 



ROW 

8 

0. 




0 * 

0 . 

0* 

0. 

0 . 

0. 

0 . 


0 . 

0. 

0* 

0. 

0. 

0. 

0 . 


100 


ooooooo ooooooo 


ROW 9 0. 

ROW 10 0. 

ROW II 0. 

ROW 12 0* 

ROW 13 0. 

ROW 14 0. 
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INPUT MATRICES FREQUENCY 2 



— 

• — MATRIX 

M A 

DIMENSIONED < 

4X 6) 

ROW 

1 

0. 

— .3000E + 01 

. 7000E+01 

• 1000E+01 

ROW 

2 

0. 

.60006+01 

.90006+01 

.10006+01 

ROW 

3 

0. 

.50006+01 

.30006+01 

-.10006+01 

ROW 

4 

0. 

.2 0 006 + 01 

1000E+01 

•3000E+01 

— 

— 

MATRIX 

M5 1 

DIMENSIONED I 

4X 61 

ROW 

1 

. 30CCE+C0 

•60006+01 

.60006+01 

.40006+01 

ROW 

2 

-.6000E+00 

.50006+01 

• 3000E+01 

•3000E+01 

ROW 

3 

-.5000E+00 

.20006+01 

. 1000E + 01 

•2000E+01 

ROW 

A 

-.2000E+00 

.20006 + 01 

0. 

.20006+01 

— 

— 

MATRIX 

C3 I 

DIMENSIONED ( 

8 X 1) 

ROW 

1 

•3000E+00 




ROW 

2 

. 1000E+00 




ROW 

3 

-.60006+00 




ROW 

4 

• 20006 + 00 




ROW 

5 

-.50006+00 




ROW 

6 

.10006+00 




ROW 

7 

-.20006+00 




ROW 

8 

• 3000 E + 00 




OUT 

P 

U T M A T R 

ICES FREQUENCY 2 



ROW 

l 

MATRIX 

0 . 

M4 1 

.34976+02 

DIMENSIONED 1 
. 7000E + 01 

7X 7 ) 

■2001E+02 

0 . 

0 . 

0 . 

ROW 

2 

0 . 

-.49996+02 

• 9Q00E + 01 

-♦3000E+02 

0 . 

0. 

0 . 

ROW 

3 

0 . 

.50006+01 

.30006+01 

-. 1000E+01 

0 . 

0. 

0 . 

ROW 

4 

0 . 

•20C06+01 

-. 10006 + 01 

. 3000E *01 

0 . 

0 . 

0. 

ROW 

5 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

ROW 

6 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0. 

ROW 

7 

0 . 

0 . 

0 . 

0 . 

0. 

0 . 

0 * 

ROW 

1 

MATRIX 

.35016+01 

M5 DIMENSIONED I 

•5013E+00 • 6000E+01 

7X 71 

•4000E+01 

0. 

0 . 

0 . 

ROW 

2 

-.50006+01 

-.17506+01 

• 3000E+ 01 

.30006+01 

0 . 

0 . 

0 . 

ROW 

3 

-.50006*00 

•2000E+01 

. 1000E + 01 

.20006+01 

0 . 

0 . 

0 . 

ROW 

4 

-.20006+00 

.2000E+01 

0 . 

.20006+01 

0 . 

0 . 

0 . 

ROW 

5 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

ROW 

6 

0 . 

0 . 

0 . 

0 . 

0. 

0 . 

0 . 

ROW 

7 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 



— 

MATR IX 

C 3 1 

DIMENSIONED I 

14 x n 





ROW I -.35016+01 
ROW 2 .lOOOE+OO 
ROW 3 • 5000E+0I 
ROW A .2000E+00 
ROW 5 -.50006+00 
ROW 6 • 1000E *00 
ROW 7 -.2000E+00 
ROW 8 .30006+00 
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ROW 9 0. 

ROW 10 0. 

ROW 11 0. 

ROW 12 0. 

ROW 13 0. 

ROW 14 0. 
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*•4 

»«******«< 

PROGRAM 

L219A2 

VERSION JUNE 29*77 IS FINISHED. 

♦ ♦ 


DATE OF 

RUN IS 

77/11/09. 


*♦4 

TIME OF 

RUN IS 

09.46.26. 

#* 
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APPENDIX A 


RELATIONSHIP BETWEEN INERTIA AND BODY-FIXED AXES 
FOR A STRAIGHT AND LEVEL REFERENCE CONDITION 


This appendix describes the transformation of equations of motion and load equations 
from inertia to body-fixed axis coordinates. The transformation developed is nonlinear 
(as shown in appendices B and C) but for straight and level flight, the equations reduce 
to a linear set. It is this linear transformation that is included in EQMOD. 

The differences between inertia and body-fixed axes for small perturbations about a 
straight and level reference condition is illustrated in figures 13 and 14. In the inertia 
axes (which are fixed in space), the motion of the body relative to these fixed axes is 
described by the velocity components in the direction of the fixed axes. Thus the 
velocities of the body are x', z' and O' for symmetric and y', 0' and 0' for antisymmetric 
motions, and the accelerations are x', z', O' and y', 0' and 0', respectively. In the case of 
body-fixed axes, the motion is described by the components of the velocity relative to 
the fixed inertial axes, but in the direction of the moving axes. For this case, the 
velocities are u, w, and q for symmetric and v, p, and r for antisymmetric motions. 
Because the axes are rotating, the expressions for acceleration contain products of 
linear and rotational velocities; they are u + Wjq, w - Ujq, q, and (v - Wjp + Uif), p, r 
for symmetric and antisymmetric motions, respectively. The factors Uj and Wj are the 
reference (in this case, lg) values of velocity in the x' and z' directions. Reference 9 
contains a development of the rigid body equation of motion in body-fixed axes. 

The symmetric and antisymmetric angles of attack a and /3 are given by w/Ui and v/Uj 
in body-fixed axes and z'/Uj +0 and y'/Ui + £*i0' - 0' in inertia axes. A point not 
illustrated in figures 13 and 14, is that the force due to gravity is constant in both 
magnitude and direction in inertia axes, but in body-fixed axes the weight force 
develops components along the axes as the axes rotate. 

Appendices B and C contain a discussion of the relationship between inertia and 
body-fixed axes based upon the formulation contained in reference 10. It is shown that 
for a general reference condition, the symmetric and antisymmetric equations of motion 
are coupled, and the relationship between inertia and body-fixed axis motions is 
nonlinear. For the special case of a 1-g level reference condition, the symmetric and 
antisymmetric equations are uncoupled and the inertia and body-fixed axis motions are 
related by the following: 

u = x'-Wj d’ " 

w = z'+Uj O’ 

q = 8' (Al) 

v = y'+Wj 0' -U] 0' 

P = 0' 
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EQUATIONS OF MOTION 


The rigid body equations of motion for symmetric and antisymmetric motion, assuming 
a i is small, is given in equations (A3) to (A6). The inertia axis equations are simple 
Newtonian equations, and the body-fixed versions are similar to those derived in 
reference 9 (eq. (4.14,12) and (4.14,13)). However, equations (A3) to (A6) do not assume 
that the origin of the axis system is at the c.g., nor that a j is zero (stability axes), but 
does assume a straight and level reference condition (8 o in ref. 9 is zero). These 
equations follow the practice of reference 10 in using a prime for inertia axis quantities. 

SYMMETRIC EQUATIONS - Rigid Body Equations of Motion 

Inertia Axes 

Mx + MAz 0' = F x ' 

Mz' + MAx cg d' = F z ' 

lyyO' + MAz cg x*' - MAx cg z*' + MgAz cg d' = M y » (A3) 

Body-Fixed Axes 

Mu + MAz cg q + MWj q + Mg<9' = F x 

Mw - MAx cg q - MUj q = F z 

I yy q + MAz C g u - MAx C g w + MAz cg W j q 

+ MAx cg Uj q + MgAz cg d ' = M y (A4) 

ANTISYMMETRIC EQUATIONS 

Inertia Axes 

My' + MAz cg <p’ + MAx cg \j) - F y » 

I xx 0' “ ! X y y' ' MAz cg y' + MgAz cg 0' = M x ' 

■zz p - *xz + MAx c y' = M z . 

Body-Fixed Axes 

Mv - MAz cg p + MAx cg r - MW j p + MU j r - Mg 0' = F y 

*xxP - l xz' r - MAz cg ^ + MAz cg W I P - MAz cg U l r + M g Az cg^ f = M x 

I zz ? - I xz p + MAx cg v - MAx cg Wj p + MAx cg Uj r - MgAx cg 0' = M z 


(A5) 


(A6) 
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where: 


M = mass 

I = inertia 


Ax c .g= x REF- x c.g. 
Az c. gr z R E F ' z c. g. 


AERODYNAMIC FORCES 

Aerodynamic forces are normally quoted in terms of lift, drag, side force, and rolling, 
pitching, and yawing moments (L, D, Y, 1, m, n), these quantities being defined along, 
and normal to, the relative airflow. The aerodynamic forces required in the equations of 
motion are the components of those forces and moments in the direction of the axes. 

Since the inertia and body- fixed axes make different angles to the airflow, the 
expressions for the forces and moments are different. Figure 15 shows these forces and 
moments for symmetric perturbations assuming that ax is small; that is, cos aj = 1 and 
sin eq = eq. In inertia axes, the angle between the axes and the flight path is a| 4- 
z'/U|, and in body-fixed axes it is ax + w/Ux; since w = z' + Ux#', the difference is 
clearly O' , and the following relationship exists between the inertia axis forces (X', Z[ 

My') and the body axis forces (X, Z, My): 

X = X' + Ld' 

Z = Z-D d' 

M y = My' (A7) 


Figure 16 shows the forces for antisymmetric perturbations, and by similar reasoning to 
the symmetric case: 


Y = Y-L0' 

M x = M x ' 

M z = M z '- m0' (A8) 

In the transformation of the roll and yaw moments, the small-angle approximation for 
ax has not been used. 
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(a) Inertia Axes 



M . = m 

y 



Figure 15. — Aerodynamic Forces and Moments for Symmetric Perturbations 
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AERODYNAMIC DERIVATIVES 


Equations (A9) and (A10) show the inertia axis perturbation aerodynamic forces written 
in terms of the customary aerodynamic derivatives. The nondimensionalizing constants 
are those adopted in reference 10. The derivation of these expressions from the force 
equations in figures 15 and 16 is outlined in appendix D. 


SYMMETRIC 


F x' " Q S W 


( C LA a l C Da) Uj + ( C L 1 + C L a °T C D a ) Uj 


( C La ^ C °a) 2Uj Uj + ( C LA a l + C L a aj + C D a 




( c i^“i * c l- £ 'i + **' 


(< 


+ * C L§ 01 \ ~ Cl 


: D S ) 5 + (CLf “I - C Df) 2U-, 

F z’ = ? s w[(- C D*“l " C LA^ + (- C D 0 «l " C L^D7 

+ (' c Di"i “ c l£ jiff, uj + ( c D(i“i 2 -c Da 

/ \ -2 


(- C d£“1 - c d^“i - c 4- c LAj2i7j 


( 


+ rCr^ _ Oi i - Cl 


*-*8 * ~ F 5 


c z 


M y ' qS w c C m A ^ + C m ^ ^ + C m g 2Ui ^ 


+ + ( C m * + ^7^ «' 

5 2^7 5 ] 


+ 2u 2 ^ + ^ m 6 ^ + 
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where: 


q = Dynamic pressure 
S w = Wing area 

c = M.A.C. 

b = Span 

5 = Control surface angle 


ANTISYMMETRIC 

Fy' = Q Syy 


c Zl + c * b y Jc 
_ yp U, + L yp W]TJJ + -yp 


+ C. 


“1)*' 


+ ( C ^“’ + Cy p) 2U7 + C y£ ~2 *' ~ C y(3 




(- c y0 + c yf)w^' + C v^ * 

2 & 1 5 ] 


+ C v ,{ + c* 

y s yg 


M x‘ - <1 S W b [( c £^ cosa l • c n j3 sin “l^|rj + ( C £^ cos “l - C n| sin “l^TJ^ + cosa, 


- C np sina, 


)“1 + ( C £p cosa l "1 - C n^ sina l “1 + CjACOsa, -C n A sina,^^ 


+ ^ C ^c°sa; sina ^—2 '<?' + C 6(J cosa, + sina,\ *' 

+ (- C £^ cos “l + c n£ sina I + CjACOsa, -C n a sina A jjjj *' 

( + Cj£cosa, - C n ^sina^ ~J2 *' + ( c £ s cosa l ~ c n 5 sin “l) 5 
+ (Cf,* cosa, - C„| cosa,^ ^ jj 



M z » = q S w b 


^ C n f3 cos “l + + (c n ^cosa l + C^sina,^ 

(c nj} cosa, + C^sinajj a, <p' + ^C^cosa, a, + C^sinaj a y 

CA sina,) 2lJ i l' + (c n *cos<*, + Cj* sina l) ^72 


+ C nA cosctj + C 


*(-s 


cosqj -Cg^sinaj^ \Js' + ^-C n ^cosa!j -Cg^sinctj + C nA cosa;j 


+ 


C Cf sina l) 2TTj + ( C n* cos “l + CjC-sina, ^ ^- 2 


( C n 


+ I C„. cosafj + Cg sina 


) 5 + ^ c n| cos “l + C £f sin “l) 2D] * j 


The DYLOFLEX system is provided with an option in EQMOD to replace the rigid body 
force elements calculated by the aerodynamic theory with terms calculated from rigid 
aerodynamic derivatives obtained from FLEXSTAB or any other source. The program 
will generate these rigid body force elements from the rigid body derivatives. Tables 2 
and 3 show these elements and their locations in the equations of motion coefficient 
matrices. The quantity aj in these figures should be understood to be the angle of 
attack of the coordinate system in 1-g flight. Since cq is calculated in an aeroelastic 
trim analysis, care must be exercised to ensure that cq is defined in the same way as in 
the dynamic analysis. Also, since the value of aj in the trim analysis is a function of 
the assumed structural constraint, an allowance should be made for the difference in 
constraint between the trim and the dynamic analysis. Tables 2 and 3 assume that the 
moment derivatives Cq and C m and C n are referred to the origin of the DYLOFLEX 
axis system. This may not be so, and EQMOD will have the facility to transfer these 
derivatives to the origin. 
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Table 2. — Formulation of the Rigid-Body Symmetric Generalized Aerodynamic Stiffness and 
Damping Matrix Elements Using Stability Derivatives 


Aerodynamic stiffness matrix 



x COL 

z COL 

fl COL 

5 COL 

x row 

0 

0 

^ S w( C 0 a + C D " a 1 C D A " a 1 ^ 

\ a R «£ u r u e 

-^Cl -a, C, +a?C L A +^0 \ 

a R a E U R U E ' 

5 

C 4„~“' °4 E ) 

z row 

0 

0 

^ ( C L a ~ a 1 C L C D a ‘" a 1 C D* ) 

\ a R u R at R u R / 

5S w( c 4 R ta ' c t>5 R ) 

d 

row 

0 

0 

q S W E ("C m C* ) 

\ a R u R / 

- q Sw c C m 


Aerodynamic damping matrix iMg) 

x COL Z COL 9 COL 5 COL 


\ 

^ Sw 

5 ^ , 

( Cq + Cq ~' Q! i ^L 

u 1 V % a E a R 

q 

( C D a + C D a + C D* + C d£ 

V q R “q E « E a E 


x row 1 

. . f °C)A l da 

U 1 V U R U E 

2 U, 

0 


'^ c l a - a iC L A ^ 

Ur U £ / 

* a 1 C L " C L, " C L 1 ) 

q e ] r 'e* 


~ a 1 C L A-a 1 C L a " a 1 C L* ~ a \ C l£ \ 
q R q £ a R a E ; 


z row 

5 Sw 

/ ClA +a i CqA \ 

^ V U R U R ) 

q S w 

— N r + “’ c ° 0r ) 

qs w c 

2Ui 

( C L A + C L^ + O] CqA + OiC d ^ \ 
\ q R a R q R a R 1 

0 

& 

row 

q SwE 

- C "* R 

q Swu c 

— — C 

U 1 "S 

2 

qS^ c 
2U 1 



( ~ ^m A ~ C m£ ^ 
\ q R a R / 

0 


x COL' z COL' ®COL' anc * ^COL are the column locations of the x, z, d, and 6 freedoms. 
{Note: There may be more than one control surface freedom.) 

The (M^] , [MgJ elements are defined in the inertial axis system. 
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Table 3. -Formulation of the Rigid-Body Antisymmetric Generalized Aerodynamic Stiffness and 
Damping Matrix Elements Using Stability Derivatives 


Aerodynamic stiffness matrix [M^J 


y COL *COL ^COL 6 COL 



Aerodynamic damping matrix [Mg) 

y COL 0 COL ^COL 5 COL 



The [M^j and [Mg] elements are defined in the inertial axis system 





TRANSFORMATION TO BODY-FIXED COORDINATES 


DYLOFLEX develops the dynamic equations of motion and load equations using inertia 
axes and provides an option to transform the coordinates to body axis coordinates by 
means of equations (Al) and (A2). The following equations show the application of this 
transformation to equations (A3) to (A6) for a straight and level reference condition. 
The equations that result are the same as the body-fixed equations except that the 
weight term in the x and y equations appears as a lift term. 

The method of implementing the transformation in EQMOD is shown in section 4.5. 

From Equations (Al): 

x' = u + Wj d' = u + Wj q 

z' = w - Uj d’ = w - Uj q 

O’ = q 

Substituting in Equations (A3) 

Mu + MAz C g q + MWj q = F x ' 

Mu - MAx C g q - MU | q = F z > 

I yy q + MAz cg u - MAx cg w + M ^z cg Wj q + MAx cg Uj q + MgAz cg Q' = M y f 

From Equations (A7) 

F x' = F x ~ ^ S W C Lj d ' 

F z' ~ F z 

My' = My 

Therefore: 

Mu + MAz C g q + MWj q = F x - q 6' 

Mw - MAx cg q - MUj q = F z 

I yy q + MAz cg u - MAx C g w + MAz cg Wjq + MAx cg Ujq 

+ MgAz cg 0' = M y 

Since for a lg Reference 

q s W C Lj = 

The Above Equations are the Same as Equations (A4) 
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From Equations (A2) 

y' = v - Wj 0' + Uj 0' = v - Wj p + Uj r 
0 = P 

# r */ _ • 

0 = r 

Substituting in Equations (A5) 

Mv - MAz cg p + MAx C g f - MW | p + MUjr = Fy' 

I xx P - I xz r - MAz cg v + MAz cg Wj p - MAz cg Uj r + MgAz cg 0' ■= M x ' 
hz' T ' I XZ P + MAx cg^ “ MAx cg W 1 P + MAx cg U l r = M z' 

From Equations (A8) 

p y' = F y + 3 S W C Lj 
^* X ' ~ ^ X 

M z' = M z + ^W c " C mi ^ 

Therefore: 

Mv - M Az cg p + MA Xcg r - MW j p + MUj r = Fy + qS^r 0 ' 

I XX p - I xz r - MAz cg v + MAz cg W| p - MAz cg U| r + MgAz cg 0' = M x 
I zz r - •xzP + MA Xcg v - MA Xcg W] p + MAx^Uj r = M z + qS w cC mi 0' 

Since: 

Q S W C Lj = 
and 

q c C m ^ - MgAx cg 

The Above are the Same as Equations (A6) 
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DISCUSSION 


It is shown that the simple transformation from inertia to body-fixed axis coordinates is 
a valid operation that should render the data more acceptable to groups accustomed to 
body-fixed axes and assist in comparing rigid body derivatives. This procedure, however, 
is valid only if the reference condition is straight and level flight. 

It should be pointed out that DYLOFLEX defines motions in a set of inertia axes with x 
defined in some convenient direction along the body. These axes are the inertia 
equivalent of what reference 9 (which uses body-fixed axes) refers to as "body” axes. 
Equally well, however, x could be defined along the reference flight path; in which case, 
the axes would be the inertia equivalent of "stability” axes. If these axes are used, then 
at\ =0. Provided aq is reasonally small, it can be (and usually is) neglected in 
symmetric dynamic analyses. However, because of the effect it has on roll/yaw coupling 
which, in turn, affects the dutch roll response, it should not be neglected in 
antisymmetric cases. 
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APPENDIX B 


RELATIONSHIP BETWEEN INERTIA AND BODY-FIXED AXES 

EQUATIONS OF MOTION 


Figure 17 shows the Laplace transform of the linearized general perturbation equations 

for an elastic airplane. The coordinates are: 

u p , v p> w p> Pp» Qp. r p Perturbation linear and angular velocities relative to a set 

of body-fixed axes with the origin at the airplane c.g. 

X P* ^p’, z p Perturbation displacements of the c.g. relative to a set of 

inertial axes. These axes are oriented to the horizontal 
through the constant Euler angles 

^P^P’^P Perturbation values of the Euler angles (rotations of the 

body- fixed axes relative to the inertia axes). 


The subscript 1 is used to denote the reference values of the above. 

U e g, U 0 a Perturbation values of symmetric and antisymmetric 

elastic coordinates chosen so that their contribution to the 
total linear and angular momentum of the airplane is zero. 
In other words, these coordinates are associated with 
free-free normal modes of vibration of the structure. 


In these equations: 

M Airplane mass. 


[ xx> l yy 


J77J 


ZZ’^XZ 


Airplane inertias. 


S 


Laplace operator. 


M S >Ds>M a ,Da,Ka The generalized structural mass, damping, and stiffness 

matrices for the symmetric and antisymmetric elastic 
coordinates. 


F x ,F y ,F z ,M x ,M y ,M z 

Qs’Qa 

'Hi 


External forces and moments. 


Time variant coefficients (functions of 6 j ,$ j ,4/ j ). 


120 



The equations assume that the airplane is symmetrical. The coordinates have been 
arranged into symmetric and antisymmetric groups. In each group, the first three 
equations are the rigid body Euler equations of motion, the next three are the auxiliary 
equations, and the last are the elastic equations. Since the inertia coordinates x p ', y p ', 
and z p ' do not appear in the equations of motion, it is usual to omit these equations. 
There is no coupling between the elastic coordinates and the rigid body coordinates 
since the elastic modes are free-free. The elastic deformations contribute to the external 
forces and moments however. 

The coefficients Ajj are shown in appendix C. These coefficients are functions of 0j, 
and ipi (the reference values of the Euler angles) and since 0j, 4>l> and if/i are 
time variant, the coefficients are time variant. It should be noted that for a special set 
of inertia axes where z' is aligned with the gravity vector, 0q = = 0 and the 

equations become identical to those shown in reference 10 (equations (6.2.29), (6.5.16) 
and (6.5.17)). The only time variant terms in the equations of motion are those 
associated with the gravity vector and it is usual to solve the equations of motion 
assuming that they are approximately constant. As can be seen from appendix C, the 
condition for these terms to be constant (if 0q and are zero) is: 9 j - 0| = 0. 


Even if the gravity terms are assumed constant, the transformation from body-fixed to 
inertia coordinates (given by the auxiliary equations) is in general nonlinear, since the 
coefficients are functions of 0j, (f>i , and tf/i . All the coefficients in figure 17 become 
time invariant if 6\ = <f>\ = = 0 . As shown in appendix C, this implies that Pj = 

Qj = Rj = 0 (that is, the reference flight condition is rectilinear. Since 0j, and 
«/» 1 are constant, the inertia axes may be oriented so that 0j = </>i = 0. Figure 18 

shows the equations for this case. 

Applying the transformation given by the auxiliary equations to the equations of 
motion results in the Lagrangian form of the equations shown in figure 19. 

If the reference condition is straight and level flight, then <f>o = Vj = 0 and tan 0q = 
Wj/Uj. If 0o is put to zero, the axes become inertia stability axes and the equations 
are shown in figure 20. The equations of motion in figure 20 are those generally used 
for loads and the symmetric and antisymmetric analyses are carried out independently. 
The forces and moments on the right-hand side are, of course, written in terms of the 
inertia coordinates. The auxiliary equations are used primarily to transform 
aerodynamic wind tunnel data, which is normally quoted in terms of body-fixed axis 
coordinates, into iner^ja coordinates. For example, figure 20 shows perturbation a and /3 
in both body-fixed and inertia coordinates. 
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Figure 17. —General Perturbation Equations 
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Figure 18. — Perturbation Equations for Rectilinear Reference Flight Conditions 




















Figure 19. — Lagrangian Perturbation Equations for a Rectilinear Reference Condition 



Figure 20. — Lagrangian Perturbation 









APPENDIX C 

EQUATIONS OF MOTION TIME VARIANT COEFFICIENTS 

The forces to gravity in the inertia axis systems are: 

F x ' = - Mgsin0 o 
F y ' Mg sin0 o cos0 o 
¥ z > = Mg cos0 o cos0 o 


Relating the inertia axis to the body axis, yields the forces due to gravity in the 
body-fixed axis system 

F x = cos0 cos \p F x ' + cos0 sini// Fy' - sin0 F z ' 

Fy = ^sin0 sin0 cosi// - cos0sin0^ F x ' 

+ ^ sin0 sin0 sin0 + cos0 cos 0 j F y ' 

+ sin0 cos0 F z ' 

F z = (cos0 sin0 cos 0 + sin0 sin0 j F x » 

+ ( cos0 sin0 sin0 - sin0 cos0 ^ Fy' 


\ 

+ COS0 COS0 F 7 ' 


F x = ^cos0j -sin0i 0 p ^ ^cos0j -sin0j 0 p ^ Mg sin0 o 

+ ^cos0j -sin0j 0 p ^ ^sin0j + cos0j 0 p ^ Mg sin0 o cos0 c 
- ^sin0j + cos0] 0 p ^ Mg cos0 o cos0 o 


-% ^U,0^p + ^u»0^p 

A u q = - Mg £sin0j cos0] sin0 o - sin0j sin0j sin0 o cos0 o 
- cos0 1 cos0 o cos0 o J 

A u -Mg |^cos0j sin0j sin0 o + cos0j cos 0j sin0 o cos0 o 
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Fy = - 

^sin0 j 

+ cos0j <f>pj ^sin0j 

+ cos0j 0pj ^cos 02 

sin0 j 0 p ) - 


^ cos 0 j - 

-sin0j 0pj ^sin0j + 

cos0 1 0p^J Mg sin0 o 


+ 

^sin0j 

+ cos0j 0 pj ^sin0 j 

+ cos02 0p^ ^sin02 + 

cos0 2 '/'p) 


^COS0| 

- sin0j 0p^ “ 

sin0j 0pj J Mg sin0 o 

cos0 o 

+ 

(sin0j 

+ cos0j 0p^ ^cos0 j 

- sin0 2 0p^ J Mg cos0 o 

cos0 o 

Fyp = A 

y e *p + 

A V,0^p + A V,0 ^P 




A v> 0 = -Mg j^-sin0j cos0j cos0j sin0 o + sin0| cos0] sin 02 sin0 o cos0 o - 
sin 0 j sin 0 j cos 0 o cos 0 o J 

A v -Mg -cos 0 | sin0j cos0| sin0 o - sin 0 ^ sin0j sin0 o + cos 0 j sin0j sin0j sin 0 o cos0 o 

L 1 


-sin 0 j cos 0 j sin 0 o cos 0 o + cos 0 j cos 0 j cos 0 o cos 0 o 


^V,0 


Fz = - 


sin 0 j sin 0 j sin 0 j sin 0 o + cos 0 l cos 0 j sin 0 o + 
sin 0 j sin 0 ^ cos 0 ^ sin 0 o cos 0 o - cos 0 ^ sin 0 j sin 0 o cos 0 o J 

^ cos 0 j -sin 0 | 0 ^sin 0 j + cos 0 j 0 p^ ^cos 0 j -sin 0 j 0 p^ + 

Mg sin0 o 


^sin0j + cos0j 0 p ^ ^sin0j + cos0j 0p|j 

+|^ cos0j - sin0j 0pj ^sin0j + cos0j 0 p^ ^sin0j + cos0j 0 p j - 
^ sin0i + cos0j. 0 p ^ cos0j -sin0^ ^p) J Mg sin0 o cos0 o 

+ ^ cos0j - sin0 j 0 p^ ^cos0 j - sin0 j 0p^ J Mg cos0 o cos0 o 
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F z,p 

A w,0 

Op + A w,^p • 

f A w ,0 





A W,0 = 

-Mg 

- COS0 j COS0 j 

cos0j sin0 o 

+ COS0J COS0 

1 sin 0i 

sin0 o 

cos0 o - 



cos0j sin0j 

cos0 o cos0 Q 





> 

II 

-Mg 

£ sin0j sin0j 

cos^j sin0 o 

- cos0 1 sin0j 

sin0 o 





- sin0j sin0j 

sin sin0 o 

COS0o " COS0 1 

COS0J 

COS0 O 

cos0 o 



- sin0j cos0 j 

COS0 Q cos0 o 

] 




> 

%- 

n 

-Mg 

cos0j sin0j 

sin0j sin0 o 

- sin0i cos 0j 

sin0 o 


“1 



+ cos0 j sin# j 

cos sin0 o 

cos0 o + sin0i sin0j 

sin0 o 

COS0 Q I 


For the condition that 

0 O = 0o = 0 


then 


A u,0 

= Mg 

COS0 1 

A U, 0 

= 0 


A v,e 

= Mg 

sin0 1 sin0 j 

A V)0 

= - Mg cos0i cos#} 

A V,0 

= 0 


A w,0 

= Mg 

cos0j sin0j 

A W,0 

= Mg 

sin0j cos^i 

A W ,0 

= 0 
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From equation 6.5-17 (ref. 10) 
x 1 - u p (^co S 6 ! -sin0j 0 p j ^cos0j -sin^j 0 p j 


+Uj 


^ cos0 j - sin0 j 0 p ^ 

^cos0j-sin0j 0pj-cos0j cos^j J 

^sin0j + cos0j 0 p ^ 

^sin0j + cos0] 6pj^cos0j -sin^j |cos0j -sin0j0 p 


^sin^j + cosip t ^pj j 

+ Vj £^sin0j + cos 0j 0 pj ^sin0j + cos0j 0 p ^ ^cos0j -sin^l 0p^- ^cos0j -sin0j 0 p j 
^sin^j + sin^j \p^j - sin0j sin0j cos^j + cos0j sin0j^ 

+ w p £ ^cos0j - sirup i ^sin0j + cos#! 0 p j ^cos^j -sin0j 0pj + ^sin0j + cos0j 0 p ^ 

^sin^j + cos^j j 

+ Wj j^costfj -sin^j 0 p ^sin0j + cos0j 0 p ^ ^cos0psin0j 0 p ^ + ^sin0j + cos0j 0 p j 

^sim^ + cos^j 0 p ^- ^cos0 j sin0j cos^j -sin0j sin0j^ 


Then 


x ~|^ A x'u u P + Vw W P + A x'0 0 p + A x'v v p + A x^ + A x> 


“ 

= Up COS 0 J COS 0 J 

_ - 


+ Wp [ cos*! sinfl 1 cos* 1 + sin*, sin* 

+ 0p£-Ui sin0j cosi^ j + Vj sin0j cos0j cos0j + Wj cos0j cos0j cos^jJ 
+ v pfsin0j sin0j cos^j -cos0j sin^jJ 


+ 0 p |Vj cos0] sinflj cosv^i + V } sin0j sin 0j - Wj sin0j sin0j cos^j + Wj cos0j sin^j 
+ * p [-Uj cos0| sin$| - Vj sin0j sin0j sintf/j -Vj cos0j cos0j - Wj cos0j sin0j sin0j 
+ Wj sin0j cos^jJ 


} 
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y' - Up - sintfj + cos0j 

+ Uj ^cosGj - sin# j tfpj^sini^] + cos0j 0 p j -cos0j sin^jj 

+ v p j^sin*] + cos0 } 0p^ |sin# { + cos<?j < 9 pj|sin 0 | + cos0j 0 p ^+^cos0j -sin0] 0 p ^ 

^COS0 1 -sin0 l 0p yj 

+ Vj ^sin0| + cos0 j 0p j l^sind y + cosGj 0pj^sin0j + cos0j 0 p j +^cos0j - sin0 j 0 p j 

^cos0| - sintf] 0pj_ s j n ^j s in0j sin0j -cos0j cosi^j^j 
+ Wp ^cos0j ~sm0j 0pj + cosOj tf p j |sin0 -j + cos0 -j 0 p ^ - ^sii\0| + cos0j 0 p j 

^cos0j - sfni^ ^ 0 p jj 

+ Wj ^cos0j ~sin0| 0p^ ^sinflj + cosO] Gp^ ^sin0^ + cos0j 0 p ^j -|sin0[ + cos0 j 0 p j 

^cos0j -sin0j ^pj-cos0j sin^j sin0j + sin0j cos0j 

y “ Ay'yUp + Ay' w Wp + Ay'fl 6^ + Ay'yVp + Ay'^<P^ + Ay'^ ^ 

y' = Upj^cosOj sin0jj 

+ v p ^sin0] sinG^ sin 0j + cos0j cos^l j 

+ sin ^l sin ^l + ^1 sin ^l cos ^l sin j + Wj cos0j cosGj sin0^j 

+■ Wp £cos0i sinGj sin0j -sin0j cos0jj 

+ 0p j^V j cos0j sinOj sin0| - Vj sin0j costfrj-Wj sin0j sin#} sin0j -Wj cos0j cos^jJ 

+ 0 p J\j] cos 0 j cos 0 j + Vj sir» 0 j sinflj cos 0 j - Vj cos 0 j sin 0 j + Wj cos 0 j sinG \ cos 0 j 
+ Wj sin0j sin0jj 
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No 


- up £sin0 j + cos0 j 0pj 

-U i [ sin0 1 + cos0 j - sinfl 

+ Vpj^(sin0j + cos^j 0p) (cos^i -sin^j 0pj| 

+ Vj |"(sin0j + cos0j 0 p ) (cos^j -sinflj 0 p)-sin 02 cosfljj 
+ Wp ^(cos 02 -sin0| 0p) (cos0j - sin 0p)j 
+ Wj ^(cos0j -sin^i 0p) (eosflj -sin0| 0p)-cos0j cosQ^J 


= Up [-sin# j J 

+ Wp [cDS0 | COS# 

+ e p ^ — u 2 cos 0 2 - V j sin 0 2 sin# j - W 1 cos 0 j sin# 2 ^ 
+ Vp ^sin 0 j cos#ij 

+ 0? [V 2 cos0j cos# j -Wj sin0j cos#jj 

+ + p [01 

= ^z’u u p * ^Z f W w p + ^z'0 * A z'v V p + ^7*0 


However, for rectilinear flight, the following must be true 


01 = Pj + (^Q\ sin0i + Rj cos02^ tan#j 

o 

0 1 — Q 1 cos 02 - Rj sin 0 i = 0 

i^l = sin 0 i + Rj cos 0 i ^secflj = 0 

o°* P j = Q 1 ~ R j — 0 
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APPENDIX D 


DERIVATION OF PERTURBATION AERODYNAMICS FORCES 
FOR THE INERTIA AXIS SYSTEM 

The following shows the symmetric perturbation forces expanded in terms of the 
derivatives of the forces X', V and My with respect to the inertia axis coordinates (Q). 
Note that 5 is a control surface rotation. The force derivatives are found in terms of the 
aerodynamic derivatives using the expressions in figure 15 and neglecting products of 
small perturbations. It should be noted that the aerodynamic derivatives are those used 
in reference 10, and equations (A-l) are used to transform the coordinates from body to 
inertia. Also, the reference drag coefficient Op is understood to contain thrust and 
consequently is zero. 


The perturbation forces are expanded in terms of the inertia axis coordinates (Q) thus 


c '_^r ax 'n- ax' •' x ax' A ax'.:, . ax' . ax' x, 

F x -L-aQ-Q'-jfv x + irr 2 + »r z + w 8 + -rr~ 8 


dx 


3z 


a T d 6 ‘ - 30- 

. ax' . ax' . ax: s 

+ ^r 6 nr 5 ir 5 


F z q ^ ax /X + Etc ' 

m' a My 

K = £ 3Q~~ Q = JiT *'+ Etc. 


where (from the equations in figure A3): 



9M y _ am 
9Q a Q 


132 



The lift, drag, and moment are expanded in terms of the body-fixed axis coordinates 
thus 


( ii yjj c W 

C Lj + C L£ I fj + C L a r \T l + C L* 2Uj TjJ 


c 2 


+ C 


L q 2Uj 4 T l; 


q + C L^ ITT 2 q + c Lx 5 + C l£ 


q 2U 


5 2U 


■.*) 


d - 


( c /? c ^ 


+ Etc. 


m qS w c + c mQ Uj + EtC ' 


or - substituting from equations (Al): 


u = x' - Wj 6 ' 
w = z' + Uj 6 ' 
q = Q' 


The expansion in terms of the inertia axis coordinates is: 

L - 5 s w [ C L, + c Lg % + % % + ^L S 2^ '%[ + (- c Lq «; + c La ) e' 

+ ( C L§ 2U| + C L§ 2Uj) ® 


7 t 2 

q 2U 2 


+ C LS +C L S « +C LJ 2^ S . 


D = <r % [ggj + C Do ^ + Etc. 

m = q S w c [C mi + C m „^ + Etc. 

. W, 

“i = U7 
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Finally 


- [fr(“i + S7) +Lx °-i] x ' + [§(“’ + %) +Lx 
+ [lf(“ + ^) +Lx °-w} z " + |f("i + u^) + Lx0 * 


3D 

w 


V, d ' , 5, and 5 


Neglecting products of small perturbation such as z' x', z' z', etc 
in the second term in the expansion of L z r : 

c .= fi£w „ q s w 1 

F * [ U, C Lq“1 U, Cd oJ 


. fi s w r .1% „ q % „ . 

+ [— C L„“1 + Uj C Lr U, C D a J Z 


q S\v c q s w c 

T~ C T ^1 - 5 — Crj* 

2Uj L~ 2U \ D c* 

q S W ( C L a - a l C L G ) a l - q S w ( C D & - a l C 
5 l^-( C L & + C L 9 )“l-q S W2U^ ( C D* 


+ Ci 


q s w c2 „ 

2U, C L 5 


<l% c2 „ 1 V 

“1" - 2Uj' C D§ 9 


q s w c l 5 “i _ q s w c d 5 8 


qS W 2U, C Lj ' q Sw 2U, C D| 8 ’ 


and similarly for F z ' and 



Jfl' + similar terms in 


and those which arise 



134 



REFERENCES 


1. Hiller, R. D. ; Kroll, R. I.; and Clemmons, R. E.: Dynamic Loads 

Analysis System (DYLOFLEX) Summary. NASA CR-2846-1, 1979. 

2. Miller, R. D. ; Richard, M. ; and Rogers, J. T. : Feasibility of 

Implementing Unsteady Aerodynamics Into the FLEXSTAB Computer 
Program System. NASA CR-132530, October 1974. 

3. Kroll, R. I.; and Clemmons, R. E. : A Computer Program To Generate 

Equations of Motion Matrices, L217 (EOM) , Volume I: Engineering 

and Usage. NASA CR-2851, 1979. 

4. Miller, R. D. ; and Anderson, L. R. ; A Program for Calculating Load 

Coefficient Matrices Utilizing the Force Summation Method, L218 
(LOADS). Volume I: Engineering and Usage. NASA CR-2853, 1979. 

5. Dusto, A. R. ; Hink, G. R. ; et al: A Method for Predicting the 

Stability Characteristics of an Elastic Airplane. NASA CR-114712 
through 114715, Volumes 1 through 4, 1974. 

Volume I FLEXSTAB Theoretical Description 

Volume 2 FLEXSTAB 1.02.00 User’s Manual 

Volume 3 FLEXSTAB 1.02.00 Program Description 

Volume 4 FLEXSTAB 1.02.00 Demonstration Cases 

and Results 

6. Shah, P. C.; and Heidergott, K. W. : Linear Systems Analysis Program, 

L224 (QR). Volume I: Engineering and Usage. NASA CR-2861, 1979. 

7. D’Auria, Patrick: DYLOFLEX Modifications to FLEXSTAB Volume II: 

User’s Manual, Volume III: Programmer’s Manual. NASA CR-2863, 1979. 

8. Miller, R. D.; and Graham, M. L. : Random Harmonic Analysis Program, 

L221 (TEV156) . Volume I: Engineering and Usage. NASA CR-2857, 1979. 

9. Etkin B. : Dynamics of Flight. John Wiley and Sons, Inc., 1959. 

10. Dusto, A. R. ; Hink, G. R. ; et al: A Method for Predicting the 
Stability Characteristics of Control Configured Vehicles. 
AFFD-TR-74-91 , November 1974. 


* U.S. GOVERNMENT PRINTING OFFICE: 1979-635-004/28 

135 


NASA CR-114712 
NASA CR-114713 
NASA CR-114714 
NASA CR-114715 




I RpfiOK ho 

NASA CR-2855 


4 Title anti Suhiitle 


2 Government Accession ho 


EQUATION MODIFYING PROGRAM, L219 (EQMOD) 
Volume I: Engineering and Usage 


3 Recipient's Catalog Nu 


5 Reoon Date 

October 1979 


6 Performing Organization Corte 


7 Aulhor(s) 

R. D. Miller, R. J. Fraser, M. Y. Hirayama, 
and R. E. Clemmons 


9 Performing Organization Name and Address 

Boeing Commercial Airplane Company 

P.O. Box 3707 

Seattle, Washington 98124 


12 Sponsoring Aqenev Name and Address 

National Aeronautics and Space Administration 
Washington, DC 20546 


IS Supplementary Notes 


8 Performmq Organization Report No 

D6-44464 


10 Work Unil No 


1 1 Contiart or Grant ho 

NASl-13938 



14 Sponsoring Agency Code 


Langley Technical Monitors: Robert C. Goetz and Boyd Perry III 

Topical Report 


16 Abstract 

L219 (EQMOD) is a digital computer program available for execution on the CDC 6600. The 
program modifies matrices according to card input instructions and prepares files of matrices 
suitable for use in the Linear Systems Analysis program (QR) and the Random Harmonic 
Analysis program. L221 (TEV156). 

The particular field of application of the program is the modification oi the theoretical equations 
of motion and load equations generated in DYLOFLEX by the Equations of Motion program 
(L217) and the Load Equation program (L218), respectively. 

Program usage and a bnei description of the analysis used are presented in volume 1 of this 
document. Volume 11 contains a description of the design and structure of the program to aid 
those persons who will maintain and/or modify the program in the future. 


17. Key Words (Suggested by Authorlil 1 

Equations of motion 
Active controls 
Load equations 
Stability derivatives 


19. Security Classd lof this report) 

Unclassified 


18. Distribution Statement 


Unclassified - Unlimited 

Subject Category 02 


22 Price" 

$7.25 


20. Security Classif. (of this page) 

21 No. of Pages 

Unclassified 

137 


‘For safe by ihe National Technical Information Service. Springfield, Virginia 22151 


















